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A great variety of systems in nature, society and technology — from the web of sexual contacts to the Inter- 
net, from the nervous system to power grids — can be modeled as graphs of vertices coupled by edges. The 
network structure, describing how the graph is wired, helps us understand, predict and optimize the behavior 
of dynamical systems. In many cases, however, the edges are not continuously active. As an example, in net- 
works of communication via email, text messages, or phone calls, edges represent sequences of instantaneous 
or practically instantaneous contacts. In some cases, edges are active for non-negligible periods of time: e.g., 
the proximity patterns of inpatients at hospitals can be represented by a graph where an edge between two in- 
dividuals is on throughout the time they are at the same ward. Like network topology, the temporal structure 
of edge activations can affect dynamics of systems interacting through the network, from disease contagion on 
the network of patients to information diffusion over an e-mail network. In this review, we present the emergent 
field of temporal networks, and discuss methods for analyzing topological and temporal structure and models 
for elucidating their relation to the behavior of dynamical systems. In the light of traditional network theory, 
one can see this framework as moving the information of when things happen from the dynamical system on the 
network, to the network itself. Since fundamental properties, such as the transitivity of edges, do not necessarily 
hold in temporal networks, many of these methods need to be quite different from those for static networks. 
The study of temporal networks is very interdisciplinary in nature. Reflecting this, even the object of study has 
many names — temporal graphs, evolving graphs, time-varying graphs, time-aggregated graphs, time-stamped 
graphs, dynamic networks, dynamic graphs, dynamical graphs, and so on. This review covers different fields 
where temporal graphs are considered, but does not attempt to unify related terminology — rather, we want to 
make papers readable across disciplines. 
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I. INTRODUCTION 

To get an overview of a large, integrated system, one needs 
to zoom out from the details. For many systems, from the 
Internet to the metabolism, from the proteome to the web of 
sexual contacts, an easy way of doing this is representing the 
system as a graph. A graph is a mathematical object consisting 
of a set of vertices, the units of the system, and a set of edges, 
the pairs of vertices that are interacting with each other. Usu- 
ally, such networks are the infrastructure of some dynamical 
system — data-packet traffic on the Internet, disease spreading 
on social networks, etc. — and this dynamical system is what 
we are really interested in. The advantage of modeling the 
system as a graph is that we can say much about the behavior 
of the dynamical system without studying the actual dynam- 
ics at all. We can estimate how much one part of the network 
influences another; how well the network is optimized with 
respect to the dynamical system; which vertices play similar 
roles in the system's operation; and so on l9l l38l l65l II 121 . 
Sometimes such a crude modeling framework can be made 
more powerful if one extends it to include additional levels 
of detail, for example edge weights in weighted networks 0, 
or the position of vertices in spatial networks IfTOll . In this 
review, we consider an additional dimension — time — and dis- 
cuss temporal networks, where the times when edges are ac- 
tive are an explicit element of the representation. Until re- 
cently, in most network studies, the time dimension has been 
projected out by aggregating the contacts between vertices to 
(sometimes weighted) edges, even in cases when detailed in- 
formation on the temporal sequences of contacts or interac- 
tions would have been available. Sometimes the solution has 
been to segment the data into adjacent time windows where 
contacts are aggregated into edges, and then study the time 
evolution of the network structure in these windows. Such 
an approach does not cover all aspects of the temporal struc- 
ture of contact patterns. For example, the edges between ver- 
tices of temporal networks need not be transitive. In static 
networks, whether directed or not, if A is directly connected 
to B and B is directly connected to C, then A is indirectly con- 
nected to C via a path over B. However, in temporal networks, 
if the edge (A,B) is active only at a later point in time than 
the edge (B,C), then A and C are disconnected, as nothing can 
propagate from A via B to C (Fig. [TJ. Thus, the time order- 
ing can matter a lot, and as we shall see below, the timings of 
connections and their correlations do have effects that go be- 
yond what can be captured by static networks. Accordingly, 
the main focus of this review is on methods that do not ignore 
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FIG. 1 : Illustration of the reachability issue and the intransitivity of 
temporal networks (more specifically a contact sequence). In (a), 
the times of the contacts between vertices A-D are indicated on the 
edges. Assume that, for example, a disease starts spreading at vertex 
A and spreads further as soon as a contact occurs. The dashed lines 
and vertices show this spreading process for four different times. The 
spreading will not continue further than what is indicated in the t = oo 
picture, i.e. D cannot get infected. However, if the spreading started 
at vertex D, the entire set of vertices would eventually be infected. 
Aggregating the edges into one static graph cannot capture this effect 
that arises from the time ordering of contacts. Panel (b) visualizes the 
same situation by showing the temporal dimension explicitly. The 
colors of the lines in (b) matches the vertex colors in (a). 

the consequences of the time ordering by e.g. projecting out 
the interaction times. 

When one studies a network, it is usually not the net- 
work itself (the vertices and edges) that is the object of study. 
Rather one wants to investigate a dynamical system on the net- 
work. In traditional network modeling one separates the un- 
derlying static network and the dynamical system on the net- 
work. Compared to this picture, temporal network approaches 
moves information about when things happen from the dy- 
namical system to the network, the underlying structure on 
which the dynamics happen. Systems suitable to be modeled 
as temporal networks are everywhere. The flow of informa- 
tion via e-mail messages, mobile telephone calls, and social 
media is one such system that has recently attracted much 
attention. Likewise, detailed understanding of the spreading 
dynamics of some electronic and biological viruses calls for 
taking the properties of the underlying contact sequences into 
account. Studies of many networks in the life sciences — from 
activation sequences of genetic regulation to time-domain fea- 
tures of functional brain networks — may benefit from the tem- 
poral graph approach. Food webs and other networks of 
species evolve in time with environmental conditions that are 
to some extent a result of which species are present. This type 
of feedback fits the temporal-network framework. Another 
example is self-assembled networks of wireless devices and 
other distributed computing systems. 

In general, when is temporal networks a suitable framework 
for analysis and modeling? lust like for static complex net- 
works, the system under study should consist of agents that in- 
teract pairwise, so that the interactions have both some degree 
of randomness and some regularity (i.e., there is some struc- 
ture). We also need to require similar properties for the tem- 
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FIG. 2: The limits of applicability of aggregated contact sequences, in the context of spreading dynamics. Panel (a) shows a schematic 
contact sequence (similarly to Fig. [T|b)) that would be fairly well modeled as an aggregated weighted graph assuming a standard random 
contact process, as seen to the right. Panel (b) displays a real-world contact sequence involving two vertices in a real network of escorts and 
sex-sellers 1 1 30 1 . The vertical lines show the times when the individuals are active in the data, while a line connecting the two individuals 
indicates a contact between them. Both the behavior of the individuals and the activity of the edge between them are bursty, with periods of 
intense activity followed by silent periods, (c) shows a hypothetical dynamics where one of the individuals in (b) gets a dose of, for example, 
a pathogen from the other individual at every contact, and the concentration of the pathogen decays exponentially. If the individual becomes 
sick when the pathogen concentration reaches a threshold (the horizontal, dashed line), then bursty dynamics would bring the level over this 
threshold. On the contrary, for more regular contact dynamics such as those in panel (a), it would have time to decay below the threshold. 



poral structures — they should not be too random or too regular 
in order to fit the framework. On one hand, one will always 
lose information when projecting a temporal network struc- 
ture to a static graph (see Fig. [T]for an illustration). On the 
other hand, in some cases, this loss of information is probably 
too insignificant to make up for the more complicated anal- 
ysis and modeling needed for the temporal graph approach. 
See Fig. |2|a) for an illustration of a contact sequence that is 
fairly well modeled by a weighted graph with the assump- 
tion that contact times are random, with a frequency propor- 
tional to the edge weight. The dynamical system of interest 
on the network matters too — different systems can respond 
differently to a specific temporal structure. For a thought ex- 
periment, consider the empirical bursty contact pattern plot- 
ted in Fig.|2|b). Assume that a contact triggers an increase of 
something (say, the concentration of a virus in the blood) in 
one of the vertices involved, which then decreases exponen- 
tially. Further, assume that the person gets sick and infectious 
if the virus concentration reaches a critical level. Then, bursty 
edge dynamics [T] [39] [76j [84] Q 63 1 could be of crucial impor- 
tance for that something to propagate through the network. In 
a situation with more uncorrelated or evenly distributed times 
of contact, the virus concentration would have time to fall be- 



low the dangerous level between the contacts. Thus for such a 
dynamical system, bursty edge activity would play a far more 
important role than for a system where the dynamics can be 
modeled as a branching process |75|, as is the case for many 
network-based models of disease spreading. 

A special case of the requirement that a system should have 
temporal structure for it to suit a temporal-network frame- 
work, relates to time scales 11431 . If the dynamical system on 
the network is too rapid compared to the dynamics of the con- 
tacts, or when edges are active, then there is no need to model 
the system as a temporal network. One example is the Internet 
where the data packets travel much faster than the topology 
changes. In summary, if the system is temporally and topo- 
logically connected in a way that affects the dynamics of in- 
terest, then temporal networks may be an optimal theoretical 
framework. 

The study of temporal networks is very much an interdisci- 
plinary field, where much of the development has been taking 
place in parallel, seemingly without much communication be- 
tween the different disciplines. This is reflected in a tremen- 
dous amount of overlapping terminology — one concept can 
easily have four or five different names in the literature. Our 
ambition is to give an overview of this research area in differ- 
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ent fields. We will not try to gather the theory into one unified 
framework. Instead, we hope that this review can help readers 
from one discipline to read and understand papers in others, 
aware of the confusing terminologies. 

Another review of contributions primarily from computer 
science can be found in Santoro et al. II 1341 and an overview 
of contributions from the network engineering community can 
be found in Kuhn and Oshman [83 1. 

In the rest of this paper, we will first discuss various real- 
world systems that can be modeled as temporal graphs. Then 
we go through theoretical developments, including measure- 
ments of temporal network structure, ways to meaningfully 
represent temporal networks as static networks, and studies of 
dynamical systems on temporal networks. 



II. TYPES OF TEMPORAL NETWORKS 

A. Person-to-person communication 

Records of electronic one-to-one communication are par- 
ticularly suitable for the temporal network approach, espe- 
cially in the context of the spreading dynamics of infor- 
mation or electronic viruses. Such data often come either 
in the form of lists of messages from one person to an- 
other at a point in time, or a dialogue between two persons 
within a time interval. The first type contains networks of 
e-mail messages ||39l |6"T1 11201 11571 . mobile phone text mes- 
sages 1 162, 169 1, and instant messages and messages in online 
forums ll58l l95l . Phone calls are not instantaneous but have 
a specific duration, and can thus be considered to be of the 
second type ll2Tll7Tl[T07llll6l [l20l. However, in many cases, 
call durations can be neglected and calls are assumed instanta- 
neous. In this context, temporal network modeling and anal- 



ysis of various temporal centrality measures (see Sect. IV 1 1 



can be used for designing strategies for containing the spread 
of malware in mobile devices [ 146 1. 



B. One-to many information dissemination 

The broadcast of information to anyone that might listen, in 
contrast to one-to-one communication, is another type of in- 
formation spreading between humans that could benefit from 
a temporal network approach. Typically, people have studied 
spreading events in blogs fl\ 1841 or microblogs (like Twit- 
ter) ll66l l87l . Liben-Nowell and Kleinberg's study of chain- 
letter e-mails concerns an intermediate form of information 
transfer, between one-to-one and one -to-many f%). These 
studies have until now focused on aggregated statistics, with- 
out much focus on the temporal effects we discuss in this re- 
view (with Ref. [96 1 as a bit of an exception, in that they also 
study response times), so further investigations are called for. 
Yasseri et al. [164| take the time dimension into account in 
an interesting way in their analysis of the circadian patterns 
of Wikipedia editorial activity: such activity patterns can be 
used to estimate the geographical distribution of editors. 



C. Physical proximity 

Proximity patterns of humans — data on who is close to 
whom at what time — are important both for understanding 
the spread of airborne pathogens and word-of-mouth spread- 
ing of information. Such temporal networks have long been 
inaccessible for large-scale studies. Rather, researchers have 
performed tedious fieldwork in some confined space like a fra- 
ternity or an office [ 160 1. Nowadays, it is fairly cheap to glean 
such information by using electronic devices. Here, the pio- 
neering study was the Reality Mining project, where students 
of Massachusetts Institute of Technology were equipped with 
cell phones whose Bluetooth devices could detect their prox- 
imity to others [37]. The SocioPatterns project has developed 
a platform that allows physical proximity measurements based 
on wearable badges equipped with radiofrequency identifica- 
tion devices (RFID) [24|; these devices have been utilized in 
measurements of dynamic and temporal proximity networks 
of patients [63 1, school children B142I . and conference atten- 
dees [140 1 . Because the human body acts as a shield for the 
proximity-sensing RF signals, such sensors only record con- 
tacts when the individuals are facing each other, and thus a 
contact can also be considered as indicative of communica- 
tion between the individuals G4l l64l 11211 11411 . Recording 
of face-to-face communication events has also been realized 
with infrared sensor devices 1 1 45 1 . Another type of large- 
scale proximity data comes from hospitals where contacts be- 
tween two patients that have been admitted to the same ward 
at the same time are recorded, sometimes including the med- 
ical staff 11991 11541 . Such data is important for studying the 
dynamics of disease outbreaks |64 141 1 (such as MRSA) in 
hospitals and also protocols for wireless ad hoc communica- 
tion 11211 . Similarly, for livestock disease, it is important to 
know the movement of animals between farms etc. ff4l 11581 . 



D. Cell biology 

There are a handful of systems in cell and microbiology that 
can be modeled as networks [119|. Not all of these are dy- 
namic enough in nature that one would benefit from modeling 
them as temporal networks. One of the systems that proba- 
bly fits the framework is the interactome — the set of molec- 
ular interactions in a cell. The vertices of the interactome 
are proteins or lighter molecules that can attach or otherwise 
connect to one another to perform biological functions. Fre- 
quently, these interactions are represented as a static graph. 
However, much of the biological functionality comes from 
the fact that the connections are not active all the time. For 
this reason, Przytycka et al. [126| believe that a "shift from 
static to dynamic network analysis is essential for further un- 
derstanding" of the interactome. There is already a body of 
literature investigating the temporal aspects of protein inter- 
action and gene-regulatory networks ll52l 1751 l92l 11511 . This 
is perhaps the most natural level for temporal network ap- 
proaches in cell biology — proteins are the workhorses of cells 
so any kind of cyclic, or otherwise dynamic, patterns have to 
be done by a shift in the interaction network. One can also 
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represent gene expression and regulatory networks as tempo- 
ral networks [91, 92 , 127, 165 1. In these, the vertices are genes 
that can be on (being transcribed) or not. Edges can be any 
of a number of functional relationships — that one gene (via 
feedback from RNAs or proteins) affect the transcription of 
another, or that two genes code for proteins that interact, or 
that they are close to each other on the DNA, etc. 

Another biological network that changes over time is the 
metabolism — the set of chemical reactions that occur in a 
healthy organism. The vertices in metabolic networks are 
molecular species that are connected if they are involved in 
the same chemical reaction. At any given time and subcel- 
lular localization, only a part of the entire biochemical reac- 
tion system is active. This situation changes with time, and 
temporal networks can potentially capture its dynamics ll26l . 
However, the change comes about via alterations in the influx 
to the cell that reflects the overall state of the organism's body, 
or is controlled by genes. Both these processes are relatively 
slow compared to the conversion of molecules, so probably 
a temporal-network analysis of metabolism does not need the 
more elaborate methods that we mention in this review. 



E. Distributed computing 

Much of the early theoretical developments on temporal 
networks come from computer science. There are many dif- 
ferent types of distributed computing systems but they all con- 
sist of fairly independent computational units spread out over 
some network fi4ll . Since the computation runs in parallel 
to the information spreading between the units, they typically 
need to operate with information that is of different age . To 
study such a system theoretically, a central problem is esti- 
mating and controlling the age of the information that is ac- 
cessible to the vertices. 



F. Infrastructural networks 

Most infrastructural networks change so slowly that there 
is no point in modeling them as temporal networks. Take the 
Internet as an example: the dynamical system in question — 
the flow of data packets — operates globally at a time scale of 
seconds. The fastest changes to the network topology come 
from new business agreements between subnetworks that are 
already in physical contact. These happen globally a few 
times per minute, but compared to the size of the entire In- 
ternet, this is so slow that one can probably safely assume 
it is static BUI 11251 . However, for some types of transport 
networks, e.g. the air-transport network, it can be meaningful 
to apply certain temporal network concepts such as temporal 
path durations and centrality |120|. In this case, edge activa- 
tion sequences correspond to scheduled transport connecting 
vertices, such as individual flights or trains. 



G. Neural and brain networks 

Networks of neural connections represent another class of 
biological networks that may benefit from the temporal net- 
work approach. There are several levels of structural and 
temporal connectivity, from the spiking patterns of individ- 
ual neurons to more coarse-grained physiological or func- 
tional connections between brain areas |20, 138 1. For the lat- 
ter, there are different experimental modalities with different 
tradeoffs regarding spatial or temporal resolution. Electroen- 
cephalography (EEG) and magnetoencephalography (MEG) 
measure electrical signals and perturbations of the extracra- 
nial magnetic fields, respectively, and have fairly good tem- 
poral but poor spatial resolutions. Functional magnetic reso- 
nance imaging (fMRI) detects changes in regional brain ac- 
tivity by measuring blood oxygenation levels, and it has a 
high spatial but poor temporal resolution. Regardless of the 
experimental technique used, the typical approach is to use 
time series associated to vertices (individual sensors for EEG 
and MEG, three dimensional regions — voxels — or larger ag- 
gregated regions for fMRI), and assign an edge between two 
vertices at a given point in time or within some time window, 
if the signals are correlated or in phase. Such networks are 
called functional networks; in this abstraction, the existence 
of a link represents simultaneous activation of brain areas, in- 
dicating a functional connection between them (see, e.g. Il20ll ). 
Naturally, such functional connections reflect the properties 
of the underlying anatomical connectivity network between 
brain areas mediated via neuronal fiber bundles (the structural 
network), but functional links may appear between brain ar- 
eas that have no (or only few) direct physical connections. In 
general, in temporal brain functional networks, the temporal 
links represent the time dynamics of simultaneous brain area 
activations - while the structural substrate network is static on 
such time scales, functional link activations vary in time. 

As an example of a static approach to functional brain net- 
works, De Vico Fallani et al. l35ll use correlations of EEG 
time series to derive a directed network in which they ana- 
lyze the motifs (over-represented subgraphs of three vertices). 
Such a study would be even richer in the temporal network 
framework where the motifs would represent temporal sub- 
networks. There are, to our knowledge, only a few papers 
where the time domain is directly taken into account: Valen- 
cia et al. II 1561 study functional brain networks reconstructed 
from MEG data with the phase-locking criterion, and show 
that the functional connectivity varies with time and frequency 
during the processing of visual stimuli, while certain network 
features such as small-world characteristics are maintained 
(see also [149]). Dimitriadis et al. [36| investigate brain dy- 
namics as measured with EEG during mental calculations, and 
identify hubs that facilitate communication in the underlying 
functional networks. Bassett et al. fl2| monitor the evolution 
of a brain network while the subject is learning a simple motor 
task. In addition, it would be of great interest to measure the 
dynamics of functional networks when the applied stimulus 
is also time-dependent, especially with naturalistic (close-to- 
real-life) paradigms such as watching a movie or listening to 
music in the fMRI scanner (see e.g. Ref. If72l ). 
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FIG. 3: A temporal network of zebras. The figure is adapted from 
Tantipathananandh et al. 11501 . The data comes from Sundaresan et 
al. 11431 . Each horizontal line corresponds to one individual. The 
contacts between individuals are not shown; instead the clusters as 
identified by the algorithm in Tantipathananandh et al. are illustrated 
by the colored squares. 




FIG. 4: Contact sequences and interval graphs. This figure illus- 
trates the two fundamental temporal network representations in our 
discussion — contact sequences (a) and interval graphs (b). The times 
of the contacts are states next to the edges. We also visualize the con- 
tacts timelines (grey bars). In these the contacts are marked by black 
bars or fields and the time lines range from t = to t = 20 (with t = 
to the left). In the former, contacts occur at points in time whereas 
the contacts are extended in time in the latter. Timeline plots like 
those in Figs.[TJb) and[2|b) are another suitable depiction of contact 
sequences, for contact graphs such illustrations are not as readable. 



H. Ecological networks 

Ecological networks capture the interactions between 
species or other categories of organisms II 1 231 11371 . They 
could be trophic, i.e. showing which species prey on another, 
or mutualistic, i.e. representing how two species engage in a 
relationship that benefits both. In many aspects, ecological 
networks are dynamic [34). As an example, they can change 
with the seasons as organisms go through different phases of 
their life cycles (and thus have different capabilities and needs 
with respect to their interaction with others) B1171 11551 . As 
another example, at longer time scales, ecological networks 
change through evolution. Since most methods of this review 
are specialized to cases where the time scale of topological 
changes of the network is not too much slower than the dy- 
namics of the network (the flux of matter between species), 
studies of evolutionary effects might not require the tempo- 
ral network approach (at least in their traditional sense, cf. 
Ref. fT66 j). For faster changes of the interaction patterns, in 
response to environmental changes, yearly and circadian cy- 
cles, etc., temporal networks could provide a useful frame- 
work. 

In population biology one also studies proximity and mo- 
bility networks of animals CH f3T] HUH [IS1 HSOl ; see Fig. [5] 
These are, just like human proximity networks, prime exam- 
ples of systems where the temporal dimension can affect dy- 
namical systems like disease and information spreading, and 
are thus apt for temporal-network analysis. In Ref. |3), dy- 
namical patterns of cattle movement were analyzed with tem- 
poral networks where vertices represent premises, and edges 
cattle movement among premises. 

I. Other systems 

The above-mentioned systems are far from the only poten- 
tial applications of temporal network modeling. Probably, 
the easiest way of finding more examples is to look at the 
complex-network literature and to ask oneself if a certain sys- 
tem has enough temporal structure for a temporal-network ap- 



proach. An early paper on time-evolving network considered 
supply networks for the manufacturing industry 1331 . There 
are likely other economic systems that would benefit from 
temporal network modeling. Networks that, like citation net- 
works [105] are normally thought of as strictly growing, could 
show temporal effects in the growth that could benefit from 
being studied in a temporal-network framework. 



III. PRELIMINARIES 

The temporal networks we consider in this review can be 
divided into two (rough and overlapping) classes correspond- 
ing to the two types of representations illustrated in Fig. |4] 
In the first representation (Fig. |4j a )) there is a set of N ver- 
tices V interacting with each other at certain times, and the 
durations of the interactions are negligible. In this case, the 
system can be represented by a contact sequence — a set of 
C contacts, triples (z, j, t) where i, j e V and t denotes time. 
Equivalently, one can represent the system by V, a set of M 
edges (pairs of vertices) E, and, for e e E, a non-empty set 
of times of contacts T e = \t\,...,t n }. Typical systems suit- 
able to be represented as a contact sequence include com- 
munication data (sets of e-mails, phone calls, text messages, 
etc.), and physical proximity data where the duration of the 
contact is less important (e.g. sexual networks). Commonly, 
authors group the contacts happening at the same discrete 
timestep into one graph (or "graphlet" in the terminology of 
Ref. ifTJI ) and present the temporal network as a time se- 
quence of graphs. Since this representation makes it tempting 
to think of the temporal-network structure as an evolving static 
network structure (which misses many of the unique points 
of temporal networks), we prefer contact sequences. Further- 
more, in many real datasets with a high time resolution there 
are only a handful of edges present at a timestep among tens 
of thousands of vertices which make the graph-sequence rep- 
resentation a little odd (but this is really just a matter of taste). 

In the second class of temporal networks we discuss, inter- 
val graphs, the edges are not active over a set of times but 
rather over a set of intervals T e = {(t\, f'j), . . . (f„, t' n )}, where 
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the parentheses indicate the periods of activity — the unprimed 
times mark the beginning of the interval and the primed quan- 
tities mark the end. The static graph with an edge between 
i and j if and only if there is a contact between i and j is 
called the (time) aggregated graph. Examples of systems that 
are natural to model as interval graphs include proximity net- 
works (where a contact can represent that two individuals have 
been close to each other for some extent of time), seasonal 
food webs where a time interval represents that one species is 
the main food source of another at some time of the year, and 
infrastructural systems like the Internet. Like for static graphs, 
it can be useful to define an index function of whether a pair 
of vertices is connected at a given time. This is the adjacency 
index (Ref. [23] call it "presence function") 

.... fl if i and /' are connected at time t 

a(i,J,t) = < n t u (!) 
10 otherwise 

Just as the largest body of literature on network theory con- 
siders simple graphs (of undirected edges that never occur 
twice between the same vertices, and never connect a ver- 
tex to itself), we put (unless otherwise stated) some further 
restrictions on the time stamps of the edges. We assume 
that a triple of a contact sequence never occurs twice, which 
means that we can order the contacts uniquely (first by the 
time stamps, then by their smallest vertex index and finally 
by their largest vertex index). Usually, we also disregard the 
order of the vertices in a contact — if contacts are considered 
directed, we will always indicate this separately. For interval 
graphs, we assume that there are no empty or overlapping in- 
tervals. More mathematically speaking, consider two intervals 
(ti, f ( '), (tj, t'.) e T e ; then the following three statements need to 
be true 

1. ti< t\ 

2. t, < ?'. 

J J 

3. t, < tj if and only if t\ < tj 

These definitions can of course be extended in many ways — 
one can think of weighted temporal networks (where a vertex 
or edge is associated with a time-dependent scalar), networks 
where the edges take some time to traverse [19| or the con- 
tacts are completed only after some duration 5t and should 
thus be represented as quadruples (/, j, t, 8t) [ 120|. Of course, 
vertices could also be active intermittently, but usually this is 
reflected in the activity of edges and we will not discuss this 
issue further. Such extensions might require modifications of 
the methods and measures we discuss in this review. Some 
such modifications are straightforward, while others are open 
research questions; several concepts of temporal graphs have 
no immediate counterpart in static graphs. 

There are other representations present in the literature. 
Those that reduce the information from the original tempo- 
ral graph by mapping them to a static graph are discussed in 
Sect. [V] A yet rarely followed path is Harary and Gupta's 
suggestion to model temporal graphs with logic program- 
ming (54). 



IV. MEASURES OF TEMPORAL-TOPOLOGICAL 
STRUCTURE 

A. Introduction 

The topological structure of static networks can be char- 
acterized by an abundance of measures (see, e.g., Il32l ). In 
essence, such measures are based on connections between 
neighboring nodes (such as the degree or clustering coeffi- 
cient), or between larger sets of nodes (such as path lengths, 
network diameter and various centrality measures). When the 
additional degree of freedom of time is included in the net- 
work picture, many of these measures need rethinking or re- 
vising. While some measures are perhaps best applied to net- 
works aggregated over chosen time periods (e.g. the time- 
dependent degree of a node can be computed as the number 
of links activated within some time window), other properties 
are directly influenced by the order of link activations. As 
an example, paths that transmit anything through the network 
need to follow time-ordered sequences of contacts, and like 
the temporal networks themselves, such paths are not static 
but change in time. In this section, we will review mea- 
sures proposed for characterizing temporal-topological struc- 
ture. Many of these build on the concept of time-respecting 
paths discussed in the first subsection, such as different cen- 
trality measures. We also address some of the few methods 
proposed for characterizing mesoscopic features and patterns 
in temporal networks; in this area, there is still a clear lack 
of methods. We conclude by discussing measures of the tem- 
poral inhomogeneities of contact sequences and information- 
theoretic aspects. 

B. Time-respecting paths and reachability 

Paths that connect nodes represent the pathways constrain- 
ing the dynamics of any process taking place on the network. 
In a static graph, a path is simply a sequence of edges such 
that one edge ends at the node where the next edge of the 
path begins (such as A to B to C to D in Fig. [TJ. In order 
for this concept to be meaningful in temporal networks, es- 
pecially in relation to dynamical processes, paths must neces- 
sarily be constrained to sequences of link activations that fol- 
low one another in time. Thus, in a temporal graph, paths are 
usually defined as sequences of contacts with non-decreasing 
times that connect sets of vertices. Kempe et al. (73 1 and other 
authors |60| call such paths "time-respecting." As an exam- 
ple, in Fig. [TJ there are time-respecting paths from A to D 
(for example (A, B, 7), (B, C, 8), (C, D, 1 1)) but none from A 
to E. In the literature, the terms "journey" |[T9l |4T1 and "non- 
decreasing path" [27 1 have also been used for time-respecting 
paths. 

The constraint of having to follow time-ordered sequences 
of contacts gives rise to differences between temporal paths 
and paths in static networks. Similarly to static directed 
networks, it might be the case that i is reachable by time- 
respecting paths from j, but j cannot be reached from i. A 
difference between directed and temporal networks is that the 
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paths are not transitive. The existence of time-respecting paths 
from i to j and j to k does not imply that there is a path from i 
to k, as seen in the example above - a path from i to k via j ex- 
ists only if the first contact on the j—k path takes place after the 
last contact on the i - j path. This is related to a fundamental 
property of time-respecting paths: they, too, are temporal, and 
begin and end at certain points in time. Thus, the existence of 
a time-respecting path that begins at i at time t' and leads to j 
does not guarantee that such a path between i and j exists for 
t > t'; in addition, a future temporal path joining i and j might 
follow a different route. Hence, the statement that "there is 
a time-respecting path between i and /' is ambiguous; such 
paths always take place within some time window. 

Thus, time-respecting paths define which vertices can be 
reached from which other vertices within some observation 
window t e [to, T]. The set of vertices that can be reached 
by time-respecting paths from vertex i is called the set of in- 
fluence of i. This is important e.g. for disease spreading, as 
it is the set of vertices that can eventually be infected if i is 
the source of infection. It may be useful to define a set of in- 
fluence at the time t as the set of vertices that can be reached 
via time-respecting paths from vertex i that begin at time t or 
later. Holme [59 1 calls the average fraction of vertices in the 
sets of influence of all vertices as the reachability ratio. 

Reversely, one can also define the source set of i as the 
set of vertices that can reach i through time-respecting paths 
within the observation window. This set consists of all ver- 
tices that can have been the source of an infection infecting 
i. Riolo et al. [ 128] points out the size of i's source set — V s 
source count — as an important quantity. Moody |108| gives 
another definition of the reachability of a vertex that increases 
in proportion to the count of time respecting paths. Again, 
as the source set is time-dependent, one may also monitor the 
source count a function of time, i.e. study how many other 
vertices may reach vertex i by time-respecting paths by time 
t', when the paths begin no earlier than t < t' . It may be use- 
ful to view the two time-dependent sets - the source set and 
the set of influence - as the past and future "light cones" for 
vertex i, i.e. the set of nodes which may have influenced i's 
current state and the set of nodes which may be influenced by 



i in the future via time-respecting paths (see also Sect. IV D 
below). 



C. Time-respecting paths with limits on waiting times 



Pan and Saramaki 1 120] note that some spreading or trans- 
port processes that follow time-respecting paths set limitations 
on the times that the paths are allowed to spend at vertices, i.e. 
times between two consecutive contacts on a path. These lim- 
itations may be from below — such that the process must be 
allowed to wait for some time before the next contact on the 
path — or from above, as for spreading dynamics, where the 
transmission has to happen quickly enough before infectious 
nodes recover, or before nodes lose their interest in forward- 
ing some piece of information. They set a limit for the latter — 
the maximum allowed waiting time at a vertex — and measure 
the reachability ratio as a function of empirical networks of 
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FIG. 5: The reachability ratio as a function of the maximum allowed 
delay of a time-respecting path at each vertex. Panel (a) shows data 
for a mobile telephone call network, (b) comes from the connections 
of an airline system. This figure is adapted from Ref. 1 120 1 . 



cell-phone calls and passenger flights (see Fig. [5]). The reach- 
ability ratio was observed to increase rather sharply around a 
characteristic time of about two days for the cell-phone data 
and 30 minutes for the airline network. These time scales can 
be interpreted as reflecting some fundamental property of the 
system. Mobile phone call sequences are known to be bursty 
and this is reflected in long inter-call intervals, and thus in- 
formation must still be further transmitted two days after its 
reception if it is to reach a large number of individuals. In- 
terestingly, 30 minutes is about the minimum allowed transfer 
time between flights for transferring passengers. 



D. Connectivity and components 

Connectivity — whether or not a pair of vertices is con- 
nected by a path — is a fundamental concept for networks. Any 
network can be divided into sets of nodes based on their con- 
nectivity; these sets, in turn, impose limitations on any dy- 
namics taking place on the network. For static networks, ver- 
tices are either connected or not, and connected components 
are defined as sets of vertices between which some path can 
always be found. As mentioned above, connectivity is not 
a symmetric relation for directed or temporal graphs. In di- 
rected graphs, the property of connectivity can be divided into 
two parts: strong connectivity, where there is a directed path 
between all pairs of vertices, and weak connectivity, where 
there is a path between all pairs of vertices if the edges are 
considered undirected. These two concepts can be general- 
ized for temporal networks. Nicosia et al. Ml 131 propose the 
following definitions: two vertices i and j of a temporal net- 
work are defined to be strongly connected if there is a di- 
rected, time-respecting path connecting i to j and vice versa, 
while they are weakly connected if there are undirected time- 
respecting paths from i to j and j to i, i.e. the directions of the 
contacts are not taken into account. On the basis of these def- 
initions, one may then define strongly or weakly connected 
components of the temporal graph as sets of vertices where 
each pair fulfills these criteria. Nicosia et al. also show that 
the problem of finding strongly connected components in tem- 
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poral graphs can be mapped into the problem of finding max- 
imal cliques in affine graphs, where an element of the ad- 
jacency matrix indicates strong connectivity between the re- 
spective vertices in the temporal graph. One should keep in 
mind that such properties always depend on the time of mea- 
surement, i.e. are in general only valid within some specified 
time window. In addition to strong and weak connectivity, 
one can define yet another type of connectivity — transitive 
connectivity — for temporal networks. A subgraph is transi- 
tively connected if time respecting paths from i to j and j to k 
implies a time respecting path from i to k. 



E. Distances, latencies, and fastest paths 

For static networks, the geodesic distance between two ver- 
tices is defined as the length of the shortest path joining them, 
path length being defined as the number of links forming a 
path. Shortest path lengths obviously influence how quickly 
anything can propagate between nodes, and their average and 
distribution determines the overall "compactness" of a net- 
work. Evidently, when the dimension of time is added to 
the picture, it is useful to define similar quantities character- 
izing how quickly vertices can reach each other through time- 
respecting paths; here, to the best of our knowledge, the earli- 
est work is by Cooke and Halsey in the 60's [ 30 1 . Here, some 
difficulties arise because the time-respecting paths are them- 
selves temporal and because the observation window is always 
finite, and choices have to be made e.g. regarding proper ways 
of averaging over quantities. In addition, the nomenclature in 
the literature has not yet converged. 

Obviously, any time-respecting path is associated with a du- 
ration, measured as the time difference between the last and 
first contacts on the path; note that some authors have called 
it the temporal path length 111201 . Analogously to the shortest 
paths that define the geodesic distance, one can find the fastest 
time-respecting path(s) between two nodes; the shortest time 
within which i can reach j is called their latency (also "tem- 
poral distance" [120|). As the concepts of temporal duration 
and link-wise distance have been used interchangeably in the 
literature, we will in the following reserve the word "distance" 
for measuring numbers of links, and "duration" and "latency" 
for measuring times (see also Sect. |IV F| below). 

The concept of latency was originally introduced in the 
study of distributed computation. A central problem in this 
area is keeping track of the age of information that a vertex 
has about other vertices. A quite reasonable assumption is 
that vertices in contact update each other's information so that 
after a contact, both vertices share the most recent informa- 
tion that either of them had before the contact. This scenario 
is similar to the SI disease spreading model that we will dis- 
cuss later, if the disease transmission probability upon contact 
is set to 100%. In terms of information spreading, this contact 
process defines the fastest possible trajectories of information 
between vertices. 

Consider the vertex i at time t in a temporal network over 
which information spreads. Then let fajij) denote the latest 
time before t such that information from j can have reached i 
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FIG. 6: The forward latency for paths from A to C of Fig. [TJ as a 
function time. There are two paths joining A and C that go through 
an arbitrary set of vertices — the first contact of the first path takes 
place at t = 7, and the duration of the path is one unit of time, i.e. 
the path arrives at C in one time unit. The next path begins at t = 11, 
and takes two units of time to traverse. If one would use periodic 
temporal boundary conditions for paths between A and C, so that the 
first observed path joining them at t = 1 repeats at t = T + 7 where 
T = 13 is the observation period limit, then the arrow would go up to 
latency 10 and then fall down linearly like for early times and thereby 
repeat a cyclic pattern. This figure is adapted from Ref. 11201 . 

by time t. We call this quantity i's view of fs information at 
time t. Furthermore, /t^fj) = t — <f>ij(J) is called fs informa- 
tion latency, or just latency, with respect to i at time t, and is 
thus a measure of how old i's information coming from j is at 
time t. Finally, the vector [</>u(l), . . . , <pi,t(N)] is called i's vec- 
tor clock. This framework was introduced by Lamport ll90ll 
and further developed by Mattern [104]. Note that the above 
definition looks backwards in time; one may also define a for- 
ward latency (called "temporal distance" in Ref. [120]) T,y(f) 
that measures how long it takes to reach j from i along the 
fastest path, when the measurement begins at time t. A bit 
reminiscent to vector clocks, Panisson et al. 1121 1 introduce 
what they call intrinsic time to denote the active time of ev- 
ery vertex. This, they argue, is useful to study the temporal 
statistics of information spreading. 

F. Average latency 

A natural use for durations and latencies is in characterizing 
the overall "velocity" of the temporal network, i.e. measuring 
how quickly vertices can on average transmit something to 
each other along the contact sequences. However, taking an 
average over the entire window of observation to get a value 
for the entire graph — or even only for a pair of vertices — is 
not that straightforward. The problem is related to the finite- 
ness of the observation period, and the fact that there are typ- 
ically different time-respecting paths between the same pair 
of vertices that begin at different points in time. Furthermore, 
close to the end of the observation window, time-respecting 
paths become rare as they do not have enough time to be com- 
pleted, and vertices may no longer reach each other. One pos- 
sible quantity for measuring the velocity of paths in general 
is to enumerate all fastest time-respecting paths between ver- 
tices and then compute the average duration of such paths. 
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However, this measure would not reflect the frequency of the 
paths, and would not be affected by waiting times before the 
first contacts of such paths: if i and k are joined by one sin- 
gle path that takes one unit of time to traverse or by ten paths 
of the same duration, this average would equal unity in both 
cases. 

For the average latency, measuring how long it on average 
takes to reach i from j, the situation is more difficult. Latency 
varies with time with a saw-tooth pattern (cf. Fig. [6]), where 
the jumps occur at points where a new fastest time-respecting 
path begins. Close to the end of the observation window, la- 
tency becomes infinite, as paths no longer have enough time 
to be completed. Averaging only over the period of finite la- 
tency is problematic: if vertices i and j were connected only 
by a single path that takes place late in the observation pe- 
riod, their average latency would be high, whereas it would 
be low if that path took place earlier. To account for this, Pan 
and Saramaki [120] proposed a pair-specific temporal bound- 
ary condition, where for every pair of vertices, the first ob- 
served path joining them is repeated once when computing 
the latency for that specific pair. Another possibility would 
be to periodically repeat the entire temporal contact sequence 
as was done e.g. in Ref. [71]. However, this procedure may 
give rise to artifacts and connect pairs of vertices that are not 
connected at all within the observation window. 

For long enough periods of observation, another difficulty 
is posed by the dynamics of vertices entering and exiting the 
system. If only edge activation sequences are observed, such 
vertex dynamics cannot generally be distinguished from edge 
dynamics — e.g. in a temporal network spanned by telephone 
calls, even if a person makes only infrequent calls, one can- 
not generally assume that the person wasn't a subscriber to 
the operator before the first call, or has left the operator after 
the last call. However, in some cases external information on 
vertices is available, and then one may choose to only include 
vertices that are known to be part of the system for the entire 
period. 

Finally, as a word of caution about the nomenclature, as al- 
ready mentioned above, some authors use the terms "distance" 
and "length" as measures of time — e.g. Kossinets et al. Il79l 
define the "distance" between two vertices as the shortest du- 
ration of any time-respecting path between them. Tang et 
al. 11471 calls the average time to reach (the reachable) ver- 
tices for time-respecting paths starting in a time window early 
in the data "temporal path length". To be fair, we should 
point out that "distance", as in the standard static graph defi- 
nition is, being a dimensionless quantity, a bit of a misnomer 
too. Furthermore, what we call average latency has several 
names: "reachability time" in Ref. [59 1, "temporal proximity" 
in Kostakos [80] 11711 . "characteristic temporal path length" 
in Tang et al. | 147 1 and "temporal distance" in Ref. 1 120 1 . 



G. Diameter, network efficiency 

There are several quantities that characterize the compact- 
ness of a static network in terms of path lengths - in addition 
to average shortest path lengths, the largest distance between 
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FIG. 7: Average latency and reachability ratio of some empirical 
contact sequences. This figure is reprinted from Ref. I59|. Each 
panel corresponds to a dataset: (a) is from contacts of an Internet 
community; (b)-(d) comes from e-mail exchange, for details see 
Holme |59|. The randomizations are Permuted Times (PT), Ran- 
dom Times (RT), Random Contacts (RC), Randomized Edges (RE) 
and All Random (AR). 



any pairs of vertices defines the diameter of the network, and 
the average over inverse path lengths of all paths the efficiency. 
One option of defining the diameter in a temporal graph would 
be to take the longest average latency (although one could 
again argue that one should avoid the mix-up between quan- 
tities whose names relate to length but that are measured in 
units of time). Again, one would then have to choose how 
to deal with infinite latencies, i.e. pairs of vertices that are 
not connected by any time-ordered path within the observa- 
tion window. Chaintreau et al. GBI gives another definition 
by requiring that the diameter should be a number as small as 
possible such that almost surely, increasing it would not make 
you find more pairs of vertices connected by a time-respecting 
path that is shorter than the diameter. The advantage with this 
definition is that it does not require all vertex pairs to be con- 
nected, which is usually the case for empirical data sets. 

Tang et al. H1471 propose network efficiency (the harmonic 
average of the latency) as a distance metric for temporal net- 
works. This measure combines the average latency and reach- 
ability ratio of Holme |59| by a harmonic mean. It was also 
mentioned in Holme |59 1 but discarded with the argument that 
the two properties (of how many vertex pairs that are con- 
nected by time-respecting paths and how fast information can 
reach between them) carry very different information about 
the function of the system and should rather not be mixed up 
(see Pig . [7j . 
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H. Minimum spanning tree 

The minimum spanning tree is an important concept related 
to paths in static weighted graphs. It is defined as a subgraph 
that is a tree of minimal total weight that connects all the ver- 
tices of the graph. Such minimum spanning trees are impor- 
tant for engineering applications as it they the cheapest way 
to reach all vertices, when the cost of an edge is measured in 
terms of weight. Gunturi et al. [49 1 presents a way to identify 
a related concept in temporal networks that they call "time- 
sub-interval minimum spanning tree". This is roughly speak- 
ing the spanning tree in an interval that has the smallest aver- 
age latency. 

As an unusual application of reachability analyses we men- 
tion Ref. |[86l that uses a temporal network approach to study 
how to navigate hot-air balloons over a planet with a pre- 
dictably changing wind field. In this case the vertices are 
cuboid volumes of the atmosphere connected if a balloon has 
the possibility of moving from one to another. (Normally, one 
probably needs more complex topologies to benefit from a 
temporal-network approach.) 



I. Centrality measures 

In network theory, numerous centrality measures have been 
defined for identifying important vertices beyond the degree, 
e.g. with respect to their average distance to other vertices 
or importance for shortest paths connecting other vertices. 
As when translating other quantities from static to temporal 
graphs, there is usually no unique candidate for the tempo- 
ral version of a static centrality measure. A rather straight- 
forward approach, however, is to replace the role of paths 
in static networks by time-respecting paths. In the words of 
Moody [ 108 1 : "time-dependent centrality measures that make 
use of the number and length of time- [respecting] paths are 
obvious choices." The closeness centrality Cc Ml 121 is for 
static networks defined as 



C c (0 = 



N- 1 

d(i, j) ' 



(2) 



where d(i, j) is the geodesic distance between i and j, i.e. the 
closeness centrality measures the inverse total distance to all 
other vertices and is high for vertices who are close to all oth- 
ers. Similarly, for temporal networks, one may be interested 
in how quickly a vertex may on average reach other vertices, 
and define the temporal closeness centrality as II 1491 



C c (i,t) = 



N-l 



(3) 



where /l,,,0') is the latency between i and j. Just like the static 
closeness centrality is undefined for disconnected graphs, 
the temporal version does not work unless there is a time- 
respecting path from j to i ending at time t the latest. This 
is, in practice, quite a hard restriction for temporal networks, 
at least if t is rather late in the observation period. In some 
situations one may want to get rid of the time dependence by 



averaging it out in Eq. Q. In this case one would have to 
choose t in the interval when is finite for all i and j, 

and integrate the sawtooth-like latency function (see Fig.[6]l in 
this interval, or to apply the boundary-condition based method 
introduced in Ref. 1 120 1 . The downside of this approach is 
that there is potentially interesting information in the intervals 
with infinite latencies, and it may therefore be better to use 
other centrality metrics. An augmented closeness centrality 
measure based on reciprocal latencies in the same vein as the 
above-mentioned "efficiency" can be defined as 



c Ei i, t) = J—Y _L_ 



(4) 



where is defined as zero if there are no time-respecting 
paths from j to i arriving at time t or earlier II 1201 1 13411 . A 
similar definition can be used for closeness centrality based 
on pairwise latencies averaged over the observation window; 
if there are no paths at all between i and j, the average latency 
is infinite for this pair and the reciprocal latency is defined as 
zero fl20l . 

The betweenness centrality Cg [38, 65] is another impor- 
tant centrality measure based on shortest paths, measuring the 
fraction of shortest paths passing through the focal vertex (or 
edge). For static networks, betweenness centrality is formally 
defined as 



C B (0 = 



Y,i*j*k Vi(j,k) 



(5) 



where v,( j, k) is the number of shortest paths between j and 
k that pass i, and v{j, k) is the total number of shortest paths 
between j and k. This definition is straightforwardly general- 
izable [ 148 ] to temporal networks by adding a dependence on 
time t and counting the fraction of shortest or fastest time- 
respecting paths that pass through the focal vertex. In the 
first alternative, the focus is on paths that contain the small- 
est number of contacts, and in the second, on paths that are 
fastest to traverse. Here, one also has to define the tempo- 
ral boundaries — options are e.g. to use the whole observation 
period, or only consider the fastest time-respecting paths that 
begin at times closest to t at each vertex j. Just like the reg- 
ular betweenness, the focus on shortest or fastest paths alone 
seems rather far from the situation in many real systems, es- 
pecially since the path durations are continuous and there may 
be paths that are only differentially slower. 

Another class of centrality measures takes its starting point 
in the assumption that something diffuses randomly around 
the network, instead of traveling from source to target along 
the shortest paths, as for closeness and betweenness. A cen- 
tral vertex in such a setup is a vertex that often is occupied by 
that something [57]. For static graphs, this approach yields 
matrix-based centrality measures like the eigenvector central- 
ity, Katz centrality and PageRank ll38lll 121 . Since this review 
focuses on methods that do not aggregate the temporal dimen- 
sion, any generalization of these measures worth mentioning 
would have to use three-dimensional tensors representing the 
temporal network. Instead of working through tensor algebra, 
we describe the algorithm of a generalization of the eigenvec- 
tor centrality. 
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1. Start with a centrality value 1 at each vertex. 

2. At every contact between vertices i and j, let the Ce 
values after the contact at timestep t be 



and 



C t +\i)=(C E (i) + (\-OC t E {j) 



CZ 1 U) = {C E U) + (l-OC E (i). 



(6) 



(7) 



When the dataset is exhausted, we have a distribution of the 
centrality that we can take as a generalization of the eigen- 
vector centrality (even though it is not derived from the solu- 
tion of an eigenvalue equation). As opposed to the shortest- 
path based measures this one becomes, automatically, time- 
aggregated. The parameter ( sets the rate of centrality trans- 
mitted at a contact. If £ is large it puts a bigger emphasis on 
recent contacts. A sensible range of £ is the interval (0, 1/2], 
but in practice it is system dependent and one needs to choose 
it by finding the value for which it ranks of the vertices as well 
as possible compared to some external data. In the upper limit, 
the vertices share their collective centrality upon a contact. 

Similarly to the eigenvector centrality applied for directed 
networks (in an iterative implementation like the one above), 
the temporal eigenvector centrality algorithm can get stuck at 
vertices. In other words, the centrality is not updated after 
the last contact, which means the vertices whose last contact 
is fairly early in the contact sequence will in the end have 
a score that has not been updated very recently. A remedy, 
similar to the PageRank or the Katz centrality, is to add a fixed 
centrality score to every vertex at every time step. To make the 
contribution uniform over time, one can normalize the score 
after every timestep (or, equivalently, increase the added value 
so that it is a constant fraction of the total centrality score). A 
version of this approach is discussed in Grindrod et al. Il47l . 



J. Persistent patterns 
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FIG. 8: Three static network quantities as a function of the sampling 
time window A (from Ref. [28)). k is the average degree (twice the 
number of edges per vertex); C is the clustering coefficient (the num- 
ber of triangles normalized to the unit interval given the number of 
connected triples of vertices, see Ref. 1 1 12 1); y is the adjacency cor- 
relation coefficient of Eq.[8] 



characteristic time scales of the behavior of the sampled peo- 
ple. This is illustrated in Fig. [8] One of the curves in this 
figure corresponds to a nifty measure of the similarity in the 
connection pattern of a vertex i by the adjacency correlation 
function (also called temporal-correlation coefficient in [ 149 1) 



7/(0 



a(i, j, t) a(i, j,t + 1) 



y/Zmut) a(i, j, t) ^Ljetm a(i, jJ+l) 



(8) 



where t represents the entire time window and (f>(i, t) is the set 
of indices j such that a(i, j, t) or a(i, j, t + 1) is one. Fig. [8] 
displays the degree k, the clustering coefficient C and y,(f) 
averaged over i and t as a function of A. 



Let us now move beyond temporal network measures re- 
lated to paths and address temporal patterns and subgraphs. 
Lahiri and Berger-Wolf Il88l discuss how to identify subgraphs 
that are persistent in temporal networks. They define the sup- 
port set S (G') of a subgraph G' (where G' Q G, is a subgraph 
of G, — the graph of all edges active at time t in a temporal 
graph) as the set of timesteps when G' c G,. They go on to 
define a "frequent" subgraph that we rather call a persistent 
subgraph as a graph that has a support larger than a certain 
threshold value. The authors have subsequently developed a 
method of detecting periodic patterns among subgraphs using 
the support set concept ||89) . 

In another study of persistent patterns in a temporal prox- 
imity network, Clauset and Eagle 1128) investigated time slices 
of an interval graph representation as a function of the dura- 
tion of the time window A. They found that the static network 
structure of the time slices depends much on A — not surpris- 
ing perhaps since the average degree grows with A — but some 
quantities shows abrupt changes of their A-scaling indicating 



K. Motifs 

Another approach for detecting significant patterns in net- 
works is related to motifs. A network motif (2) is an equiv- 
alence class of subgraphs that is overrepresented in terms of 
its cardinality with respect to some null model in a network, 
i.e. a larger number of such subgraphs can be found than in 
a randomized reference system H172L Usually, the configura- 
tion model that conserves the degree sequence but otherwise 
randomizes the network is taken as the reference system. The 
over- or underrepresentation of certain subgraphs can be re- 
lated to the function of the system, especially in directed net- 
works where the subgraphs forming motifs can be associated 
with e.g. information processing tasks. 

There are several ways to extend this concept to temporal 
networks. The easiest is to look at snapshots of the network 
taken at different points in time, or alternatively aggregated 
edges over a period of time, and count the different subgraphs 
in these snapshots. Braha and Bar- Yam [18] does exactly this 
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for an email data set, where one snapshot network represents 
contacts aggregated over a day. They look at motifs of four 
vertices and compare them to an edge rewiring null model. 
The main conclusion from Braha and Bar- Yam's study is that 
the denser motifs are overrepresented. This approach, how- 
ever, projects out the information from the order of events, 
as the motifs are no longer temporal networks themselves. 
In Ref. [5 1, Bajardi et al. apply an even simpler definition 
for their dynamical motifs that are in essence time-respecting 
paths constructed from events belonging to adjacent snap- 
shot time windows. Interestingly, when the time-reversal null 



model is applied on data on cattle movements (see Sect. VI C i 



the number of such paths is observed to be smaller than for the 
original data, indicating the presence of an "arrow of time" in 
the system. 

The communication motifs addressed in Zhao et al. 11681 
are also defined on the basis of static subgraphs, although 
temporal information is used for defining the subgraphs of 
interest. Here, the authors take an approach where commu- 
nication events in mobile telephone call records or Facebook 
wall postings are first linked to form a communication graph 
if they both share vertices and succeed one another within a 
time limit At, similarly to the path waiting-time cutoff. The 
difference to the above definitions is the absence of a snap- 
shot window; temporal adjacency of events is based on the 
time difference between two events sharing a vertex, instead 
of aggregating the entire system. Communication motifs are 
then constructed on the basis of graph equivalence of (static) 
subgraphs found in such networks. Zhao et al. observe high 
frequencies of chain, star, and "ping-pong" subgraphs that 
can hardly be found in reference networks where the event 
times have been randomly reshuffled (the Randomly Permuted 
Times null model, see Sect. | VI C\ Such motifs are seen to be 
stable over time. Zhao et al. also define "maximum flow" mo- 
tifs, where call durations are used as a filtering criterion. 

Chechnik et al. E6l define temporal motifs for gene- 
regulatory networks controlling metabolic pathways. In this 
case, the primary temporal unit is the genes, or vertices, not 
the edges. Their method proceeds by coding the time evo- 
lution of the expression level of a gene into a sequence of 
states and comparing the relative timing of the transitions be- 
tween theses states for different connected genes. Further- 
more Chechnik et al. use randomization techniques similar 
to those in Sect. |VIC| to infer statistically significant "activity 
motifs" as they call them. 

In Ref. ifSTI . Kovanen et al. define temporal motifs such 
that the equivalence classes are based on temporal subgraphs, 
i.e. the order of contacts constituting the subgraph/motif mat- 
ters. The requirements for the temporal subgraphs are based 
on the following definitions: two contact events e, and e ; are 
adjacent if they share a vertex, and they are At-adjacent if 
additionally their time difference is no larger than Af. Fur- 
ther, two events <?,- and e^ are At-connected, if there exists a 
sequence of Af-adjacent events joining e, and e^. Note that 
this sequence does not have to be a time-respecting path. A 
temporal subgraph may now be defined as a set of events 
where all pairs are Af-connected. Requiring that no events 
are skipped at nodes within the subgraph yields valid tempo- 




FIG. 9: A contact sequence (the numbers on the edges denote the 
times of contacts) involving five vertices (left), and the two maximal 
Af-connected temporal subgraphs found within this sequence when 
At = 15. After Ref. ED- 



ral subgraphs. Temporal motifs are now defined as classes of 
isomorphic valid subgraphs, where the isomorphism is taken 
to include the similarity of the temporal order of events, but 
not their exact timings. Furthermore, for every event e, there 
is a unique maximal Af-connected temporal subgraph that in- 
cludes e, and the largest possible set of events that are still 
pairwise Af-connected. Motifs based only on maximal tem- 
poral subgraphs are called maximal motifs. 

For the algorithmic solution, such temporal subgraphs may 
nevertheless be mapped into static directed graphs, whose iso- 
morphism can then be used for computing subgraph counts 
for the temporal motif equivalence classes. One can further 
merge equivalence classes by considering (potential) informa- 
tion flow within the event sequence of a temporal subgraph — 
in some cases, there are some events whose detailed order 
does not matter. When applying these definitions and a cor- 
responding algorithm to large-scale mobile call data, Kova- 
nen et al. found out that temporal motifs where the event se- 
quences may have a causal explanation are more frequent than 
sequences where the event order appears random. 

Regarding temporal motifs, there are issues that still remain 
unresolved and would warrant further consideration. A very 
important issue is that of reference or null models. For static 
networks, motif analysis is concerned with over- or under- 
representation of subgraphs when compared with a reference 
ensemble where structural correlations have been removed. 
However, what the reference ensemble should be for tempo- 
ral motifs is far from obvious, as it is evident that for tempo- 
ral graphs where the events are fairly sparse, a much smaller 
number of temporal motifs can be found in a reference en- 
semble where the event times have been randomly reshuf- 
fled. Thus any larger temporal subgraphs where the waiting 
times between events are short are almost always overrepre- 
sented, and not much information is gained from applying the 
null model. Another issue that would deserve further atten- 
tion is the recurrence of temporal subgraphs — either exactly 
the same or to some extent similar sequence of contact events 
may take place between the same set of vertices at multiple 
points in time. Understanding such recurrent mesoscale pat- 
terns might yield a lot of insight e.g. on the function of social 
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communication networks, or functional brain networks. Fur- 
thermore, for some systems, studying the time dependence of 
the occurrence counts of temporal subgraphs might also yield 
useful information. 



L. Measuring inter-contact times and burstiness 




In addition to measures characterizing the temporal net- 
works and patterns in them, it is often also useful to have 
a closer look at the smallest building blocks of such net- 
works - the vertices and edges, and the associated sequences 
of contacts. In a typical temporal network, there are multi- 
ple contacts between connected vertices, and correlations in 
the timings of such contacts play an important role e.g. re- 
garding durations of time-respecting paths and latencies be- 
tween vertices. Especially, for temporal networks of hu- 
man communication, it has been discovered that such tim- 
ings are often bursty and deviate from the more uniform 
times expected from a memoryless, random Poisson pro- 
cess 0|^l22|6D|62lin[IO2[na Overall, 
the burstiness manifests itself in broader-than-expected distri- 
butions of inter-contact times, P(r), defined either for nodes or 
their individual edges. However, because of inhomogeneities 
and broad distributions of node degree, numbers of contacts 
per node, and numbers of contacts per edge, directly interpret- 
ing the shape of P(t) can be difficult as it reflects combined 
effects of burstiness and such inhomogeneities. Because of 
this, displaying scaled inter-contact time distributions has be- 
come the norm CTO I7T1 [1071 ■ Here, nodes (edges) are binned 
according to their numbers of contacts, and the inter-contact 
time distributions are scaled by the average inter-contact time 
in each bin, P{t/t*). When compared against similar distri- 
butions for uniformly random contact times, broad tails are 
typically observed. 

Another alternative is to directly derive a measure of bursti- 
ness for any sequence of inter-contact times. Here, the sta- 
tistical properties of inter-contact times arising from a Pois- 
son process form a natural point for comparison. Goh and 
Barabasi [45 1 use as their starting point the coefficient ofvari- 
tion, defined as the ration of the standard deviation of the inter- 
contact times to their mean, cr T /m T . For a Poissonian contact 
sequence, o~ T /m T = 1. Using this quantity, the burstiness of a 
sequence is then defined as 



FIG. 10: Reachability graphs. Panel (a) shows a contact sequence 
(same as in Fig.[TJ and (b) shows its reachability graph. 



(o- T /m T -l) (cr T -m T ) 

D = = . 

{cr T lm T + 1) (cr T + m T ) 



(9) 



For finite contact sequences, the variance is always finite, and 
B € (-1, 1), such that B — 1 indicates a most bursty sequence, 
B — a sequence with Poissonian inter-contact times, and 
B — - 1 a completely periodic sequence. 



M. Entropies and other information-theoretic measures 

Information theory is the branch of science exploring the 
limits of storing, communicating and compressing data. It has 
already found applications in the network-clustering literature 



(e.g. Ref. 1 133|) and could well be extended to temporal net- 
works. One such attempt is Timo, Blackmore and Hanlen's 
entropy-based method to account for temporal uncertainties 
in communication networks 111521 . While this is not a deter- 
ministic temporal-network approach, as in most of this review, 
it points to an interesting direction for future research. Simi- 
lar ideas have been developed in the analysis of ecosystems. 
Ulanowicz define a measure "ascendency" to quantify how 
well a system process its environmental input Ml 17111551 . The 
measure, derived from information theory, can be adapted to 
ecological interactions changing in time. 



REPRESENTING TEMPORAL DATA AS A STATIC 
GRAPH 



The literature on static graphs is many times larger than that 
on temporal graphs, for a natural reason: it is usually much 
easier to analyze static graphs, especially analytically. One 
approach to analyzing temporal graphs is thus to derive static 
graphs that capture both temporal and topological properties 
of the system. The most straightforward way is to accumu- 
late the contacts over some time to form edges. This can be 
an informative approach to studying the time evolution of the 
static structure of an evolving network |22|. Typically, au- 
thors have either investigated a network-topological quantity 
as a function of time when every contact between a pair of 
vertices that has not been in contact before adds an edge to 
the graph (e.g. Ref. [58]), or they have divided the time into 
segments and studied the structure of contacts accumulated 
over those segments 11331 . As mentioned above, this trivial 
way of projecting out the temporal dimension can discard in- 
formation and this review focuses on methods that attempt to 
capture it. Nevertheless, it is useful whenever the topological 
aspects are more important than the temporal. It may also be 
possible to combine these aspects: Miritello et al. [ 107 1 have 
proposed a way of mapping the dynamic SIR model to a static 
edge percolation model, where the edge weights of a network 
whose structure equals that of the aggregated graph are de- 
fined in a way that takes into account the temporal correlations 



and inhomogeneities of edges (see Section VII I. There are 



other ways of encoding the temporal network structure into 
a static graph different than the aggregated graph that incor- 
porate more temporal information than a plain projection of a 
contact sequence or interval graph into the graph dimension. 
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A. Reachability graphs 

An alternative graph representation that might be useful at 
least for very sparsely connected contact structures is reacha- 
bility graphs, or "path graphs" [108], or "associated influence 
digraph" l27l . In such a case one puts a directed edge from 
vertices A to B if there is a time-respecting path from A to B 
(see Fig. 10 1. Such a graph, thus, shows which vertices can 



possibly affect which others. The average degree A: of a reach- 
ability graph is thus the average worst-case outbreak size mi- 
nus one. In other words, for any contact structure that supports 
a pandemic, the reachability graph will be dense (k ~ N). This 
is a drawback for reachability graphs since most methods to 
analyze networks are developed for sparse graphs. Neverthe- 
less, one can imagine very sparse connection structures that 
can benefit from such a representation. 

One work using reachability graphs is Bearman et al. Ifl4l 
that studies dating between high-school students. Their data 
fits the interval-graph framework (Fig.|£|)) where contacts can 
happen at any time during the intervals, and the only thing 
one knows is that the contact happens at least once. In such a 
case the reachability graph is not unique and the authors plot 
a maximal and minimal reachability graph. Going in the op- 
posite direction, one can prove that given a reachability graph, 
one can always construct a temporal graph that has the given 
reachability structure 11271 . 



B. Line graphs 

Although it is a bit outside the scope of this review, we 
mention that line graphs have been used to study disease 
spreading in temporal networks represented by aggregated 
time slices |97|. A line graph of a simple static graph G is 
a graph whose vertices are the edges of G that are connected 
if they share a vertex in G. Sometimes the line graph is called 
"interchange graph" or "dual graph" (the latter being a bit of 
a misnomer since the dual of the dual of G is not G). The 
point of studying line graphs in epidemiology is that they are 
closely related to the structure of concurrent partnerships — 
e.g. the number of edges of the line graph is the number of 
concurrent partnerships in the original graph. 



C. Transmission graphs 

Riolo et al. Ml 281 present a version of line graphs adding 
some more temporal information. What they call a transmis- 
sion graph thus goes beyond modeling the contact structure 
and also incorporates the disease dynamics via a parameter 6 
that quantifies the combined incubation time and duration of 
a disease. Their starting point is a graph where an edge e is 
active over an interval [f s tart(<?), f s top(e)]. There is a directed 
edge from e to <?' in the transmission graph if e and e' share 
a vertex, f star t(<?) < f s tart(e') + 5 and t^ n {e') < f sta rt(e). See the 



disease. The advantage over a pure line graph is that transmis- 
sion graphs encode the directionality arising from the order of 
non-concurrent relationships. However, in common to line 
graphs of concurrent relationships, the transmission graphs 
cannot handle edges where one vertex manages to not catch 
the disease. Furthermore, paths in transmission graphs do not 
have to be time respecting. Riolo et al. [ 128] define some dif- 
ferent versions of their transmission graphs cutting the contact 
sequence in various ways, but essentially the idea is outlined 
above. 



VI. MODELS OF TEMPORAL NETWORKS 

In general, models of networks can serve many purposes. 
They may explain the emergence of salient network charac- 
teristics, or serve to produce synthetic networks with desired, 
tunable characteristics that can then be used in computational 
experiments of dynamic processes. Another class of models is 
that of randomized reference networks, where empirical net- 
works are taken as inputs and randomization procedures are 
used to remove some of their characteristic correlations. For 
temporal networks, the number of models proposed in the lit- 
erature is still fairly limited. Below, we will first discuss some 
models proposed for temporal social networks and epidemi- 
ological contact networks, and then move on to randomized 
reference models. 

Note that the word "model" has a slightly different meaning 
when used in the context of statistical inference - here, tech- 
niques of inference and machine learning are used to extract 
statistical models, typically from limited amounts of data. 
Methods for statistical inference of temporal networks are still 
rare, although this problem is an important one for systems 
biology, e.g. in the context of temporal networks of gene reg- 
ulation, where Ref. [92 1 use a Markov Chain Monte Carlo 
approach to infer the parameters of a so-called Bayesian Net- 
work representation of gene regulatory networks [91]. This 
approach could probably have a wider use for systems also 
outside of biology. 



A. Models for temporal social networks 



I. Temporal exponential random graphs 



illustration in Fig. 1 1 The idea with the extra 5 is to identify 
the last possible time a member of an edge can transmit the 



The method of exponential random graphs (see, e.g. II 1291 ~) 
is commonly used by social scientists. The parameters of such 
graphs, inferred from empirical data, contain information on 
the importance and frequency of chosen topological elements 
and subgraphs, such as triangles and stars. A similar modeling 
framework for temporal exponential random graphs has been 
put forward in e.g. Refs. [50 53]|77]. The basic problem in 
this literature is that given an observed time evolution of a set 
of states of the vertices, representing some dynamical system 
on the graph, one should determine the parameters of a time- 
varying exponential graph model that essentially gives the 
probability of a contact to happen between a pair of vertices 
at a given time. Just like its static counterpart — exponential 
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(a) (b) time (C) 

FIG. 1 1 : Transmission graphs. Panel (a) shows an interval graph representation of a temporal network (where each edge has only one interval, 
as required by the definition by Riolo et al. 1 128 1). Panel (b) gives slightly reduced picture of the system, where a line corresponds to an edge 
in (a) and the active interval is indicated. The derived transmission graph is illustrated in (c). 



random graphs — has a connection to the Ising model, the tem- 
poral exponential random graphs can also be boiled down to 
finding a partition function, in this case a time-varying one. 
When this is done, the temporal exponential random graph 
model can be used both as a reference model for measuring 
biases in quantities of topological and temporal structure, as 
well as a generative model for tuning the structure of con- 
tact sequences for simulations of dynamical systems on top of 
these. 



2. Models of social group dynamics 

Zhao, Stehle, Bianconi and Barrat have presented a frame- 
work for modeling social networks H 1391 11671 where edges 
represent ephemeral social ties, such as being in face-to-face 
contact. Their approach is based on master (or "rate") equa- 
tions that represent the expected change in the number of peo- 
ple in a group of a certain size, and can capture observations 
like "the longer an agent interacts with a group, the less it is 
likely to leave the group; the more the agent is isolated the less 
likely it is to interact with a group." In essence, this framework 
deals with interval graphs rather than contact sequences, much 
like the neighborhood exchange model discussed below. 

In contrary, the model by Jo et al. [68 1 explicitly accounts 
for the contact timings, combining their short time scale with 
the longer time scale of network evolution. Jo et al. pro- 
pose a social network model that combines features of earlier 
social network models (focal and cyclic closure, tie strength 
reinforcement l8~5l ) with triad interactions and social interac- 
tion task execution from a priority queue, along the lines of 
Ref. Q- This model gives rise to networks with communities 
of strong ties connected by weak links; the contacts mediated 
by the links are bursty. 

B. Contact network models 

As a simple way of extending static graphs to involve a 
turnover of neighbors is was proposed by Volz and Mey- 
ers 11591 . This model is rather similar to the RE randomiza- 
tion procedure discussed below, but with the purpose of mim- 
icking the change of partnerships to generate contact struc- 



tures for disease simulations. It works by selecting two edges 
with some probability every time step and swapping them as 
in step 3 of the RE procedure. As the model is based on 
rewiring an existing network, the topology of the accumulated 
network and the contact patterns across edges have to be de- 
termined using some other models. 

In order to generate synthetic sexual contact networks, 
Kretzschmar et al. |82| have proposed a model that gener- 
ates interval graphs and allows tuning the assortativity of the 
resulting networks.. The model rules are as follows: 

1. (a) A new partnership is formed with probability p. 

The individuals participating in the pair are cho- 
sen according to the mixing function <p: 

i. draw two random individuals i and j; 

ii. decide whether they form a pair according to 

<P(i, j); 

iii. if yes, done; else go to ;. 
(b) Repeat step[Ta]A72 - P times. 

2. In every pair consisting of a susceptible and an infected, 
the disease is transmitted with probability rj. 

3. Every pair splits up with probability <x. 

The mixing function is introduced to be able to tune the mix- 
ing by degree 111 121 . It could be 

0(i,/) = W + ?pr L do) 

^max 

for assortative mixing (where there is a tendency for high- 
degree vertices to connect to other high-degree vertices and 
low-degree vertices to low-degree vertices), or 

(jfcj - k:) 2 
^raax 

to create disassortative mixing where high-degree vertices 
tend to connect to low-degree vertices. The parameter £ sets 
the strength of the assortativity or disassortativity. k m . dx is an 
upper limit of the degree — this was before the finding that sex- 
ual networks have power-law degree distributions (Liljeros et 
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al. [98 1). A modern approach would probably be to draw de- 
grees from a distribution and then connect them with a func- 
tion like <p. Except controlling for assortativity, one can also 
choose <p to create serial monogamous relationships by letting 
= 1 if ki = kj = 0, otherwise <p = 0. In a paper on measur- 
ing concurrency (the temporal overlap of partnerships) Morris 
and Kretschzmar [ 109 1 used the parameter values: N = 2000, 
p = 0.01 /day, cr = 0.001/day and 77 = 0.1 /day. We note 
that, as a model for sexual networks this model fails to cap- 
ture that overall sexual activity decreases with the number of 
partners 111 141 . 

A fair amount of subsequent works has followed ideas sim- 
ilar to the stochastic pair-formation model. We will not review 
all of them but mention the "dynamic random graph models 
with memory" of Turova [ 153 1, which effectively models the 
same type of time-evolving graphs as above, but use a frame- 
work more tractable for analytic calculations. Other recent 
modeling evolving interval graphs but including more struc- 
ture are presented in Refs. H135M136tt . 



C. Randomized reference models 

For static networks, a common way of assessing the impor- 
tance, unexpectedness, or over/underrepresentation of topo- 
logical features of empirical networks is to compare the fea- 
tures against some reference model where the network is ran- 
domized. The most widely applied reference model is the con- 
figuration model, where the links of the original network are 
randomly rewired pairwise. This reference model preserves 
the original degree sequence but yields otherwise maximally 
random networks. Then, one can assess the significance of 
quantitative topological characteristics of the empirical graph 
by either a direct comparison to averages in the randomized 
reference ensemble, by computing Z-scores for the character- 
istics with respect to their distributions in the reference ensem- 
ble given by the reference model, or by measuring the extent 
to which the dynamics of some processes differ when run on 
the empirical networks and the reference ensemble. 

For temporal graphs, a similar approach can be applied: in 
this case, the original event sequences are randomized or ran- 
domly reshuffled in a chosen fashion to remove time-domain 
structure and correlations. However, there are several kinds of 
possible temporal correlations and several time scales where 
the correlations are important, and thus no single, general- 
purpose null model can be designed (the temporal config- 
uration model). Rather, by designing appropriate null mod- 
els, one may switch off selected types of correlations in order 
to understand their contribution to the observed time-domain 
characteristics of the empirical temporal network. Such tem- 
poral null models have also been applied for studying the ef- 
fects of various kinds of correlations on dynamical processes 
(such as spreading) on temporal graphs. A typical use for 
such models in studies of dynamical processes would be to 
essentially apply all of them, and by monitoring how the dy- 
namics of the process depends on the reference models, to 
pinpoint the role of different temporal and topological corre- 
lations on the process - if removing a certain type of correla- 



tions changes the dynamics a lot, then obviously those play an 
important role for the dynamics. 

Below, we review temporal null models introduced in the 
literature, essentially following Holme |59| and Karsai et 
al. ff\\ . For this section, we assume that the temporal net- 
work is a contact sequence. Some of the methods work for 
interval graphs too; others can be modified to interval graphs 
quite straightforwardly. In the end of the section, we sum- 
marize and provide some guidelines for choosing reference 
models. 



1. Randomized edges (RE) 

This method is similar to the configuration model for static 
graphs mentioned above, with the additional ingredient that 
contact sequences of edges follow the edges when these are 
rewired. Algorithmically, the method is defined as follows: 

1 . Go over all edges sequentially. 

2. For every edge (2, j), pick another edge (?, /). 

3. With a probability 1/2 replace (/, j) and (/',/) by (2, /) 
and otherwise replace them by (2, /') and (J, /). 

4. If the move in step [3] created a self-edge or multiple 
edge, then undo it and start over from step [2] 

The times of contact over an edge are kept constant. Note that 
the two alternatives in step [3] where one is randomly selected 
are needed to remove spurious correlations if the data struc- 
ture that is used returns the vertices of an edge in a specific 
order; otherwise one would keep the number of times a vertex 
appears in the first argument conserved, which in practice can 
give quite big differences for empirical graphs, whether the 
graph is small or not. To speed up the process, one can skip 
edges that already have been rewired in step [T] (by being se- 
lected in a previous step [2}. On the other hand, this procedure 
is linear in M and rarely a computational bottleneck. 

This null model can be used to study the effect of the net- 
work topology, that is, the wiring diagram of the original net- 
work. The model also assumes that it is the edges rather than 
the vertices that govern the times of contacts — after the ran- 
domization procedure, both the numbers and timings of con- 
tacts for each vertex will have changed; however their degrees 
in the aggregated network are retained. As the contact se- 
quences follow their edges when rewiring, all temporal corre- 
lations and inhomogeneities associated with individual edges, 
such as burstiness and the distribution of inter-contact times of 
edges, are retained, as is the overall event rate at every point 



in time. The RE procedure is illustrated in Fig. 12 



2. Randomly permuted times (RP) 

As a temporal counterpart to the configuration model, one 
can permute the contact times randomly while keeping the 
network structure and the numbers of contacts between all 
pairs of vertices fixed. Technically, this is much simpler than 
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FIG. 12: Illustration of two types of randomization null-models for contact sequences, (a) shows a contact sequence (the same as in Fig. |TJ. 
In (b) it is randomized by the Randomly Permuted times procedure such that contacts happen the same number of time per edge, and the 
aggregated network topology is the same. In (c) the contact sequence in (a) is randomized by the Randomized edges (RE) procedure. With 
RE, the time sequence of the contacts along an edge is conserved, and so is the degree sequence of the original network, but all other structure 
of the topology is destroyed. (The latter statement is perhaps not so well illustrated by this figure as there are not so many graphs with the 
degree sequence of the original, aggregate graph.) 



applying the edge rewiring scheme discussed above and it 
only requires randomly exchanging the time stamps of all con- 
tacts, or just randomly reshuffling the order of time stamps in 
an array or vector. No checks similar to step 4 of the RE 
rule need to be performed for contact sequences; for interval 
graphs, one should check for non-overlap. As this null model 
retains all network structure and the number of contacts for 
each edge, its application can be used to study the effects of 
the detailed order of events, including burstiness, inter-contact 
time distributions of contact sequences on vertices and edges, 
and correlations and triggering effects between contacts on ad- 
jacent edges. The model also retains the overall rate of events 
in the network at every point in time, such as daily or weekly 
patterns in communication networks. An illustration of RP 



can be found in Fig. 12 



3. Randomized edges with randomly permuted times (RE + RP) 

This null model works as follows: first, the network struc- 
ture is randomized using the RE procedure. Then, the time 
stamps of all contacts are reshuffled with the RP scheme. Thus 
the outcome is a temporal graph, where all structural correla- 



tions and all temporal correlations with the exception of the 
overall rate of contacts (such as daily/weekly pattern) have 
been destroyed. 



4. Random times (RT) 

The RP ensemble conserves the set of times of the origi- 
nal contact sequence. Hence, although it destroys burstiness 
of events on individual vertices and edges as well as correla- 
tions between events such as triggered chains, the aggregated 
rate of events in the network is unchanged and will still fol- 
low the typical circadian and weekly patterns of human activ- 
ity EH ED EE] HSU HSU ESI- So far, results indicate ED 
that such overall modulation of the event rate does not appear 
to matter at least for the simple SI spreading model, whose 
dynamics is dominated by edge burstiness (see Jo et al. Il67l 
for a method that removes the contribution of such circadian 
patterns). However, there are cases where such patterns might 
play a role, such as more complicated spreading models (e.g. 
SIS or SIR). The random times (RT) null model destroys such 
patterns: in this approach, each contact on each edge is as- 
signed a random time within the observation time window of 
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the original data, and thus the network structure and the total 
number of events on each edge are conserved. When compar- 
ing the results of any dynamics taking place on the RP and 
RT models, the difference should then be due to the network- 
level temporal patterns of the empirical data. Note that an- 
other alternative for uniformly random contact times is gen- 
erating them from any chosen distribution or process, such as 
the Poisson process [71 1, with parameters set up so that the 
numbers of contact per each edge are on average conserved. 

The above approach assumes the degree distribution is ex- 
ternal, as it is only investigating ensembles where the topol- 
ogy and the degree of each vertex is conserved. Another op- 
tion would be to run the dynamics on Poisson random graphs 
(Newman 1 1 12 1) with the same number of vertices and edges 
as the original graph and time stamps distributed as the origi- 
nal data, permuted like RP and random like RT. 

5. Randomized contacts (RC) 

Here, one keeps the graph topology fixed but redistributes 
the contacts randomly among the edges. After this random- 
ization, the number of contacts per edge follows the binomial 
distribution rather than some broad, right-skewed distribution 
as is typical for empirical data. This randomization can be 
utilized in testing the effect of the distribution of the number 
of contacts per edge in combination with the order of events. 
Hypothetically, one would like to test the effect of the distri- 
bution of the number of contacts alone, keeping the structure 
of the temporal order of the real data. For example, a vertex 
that is active primarily in the early stage of the data would be 
so in the randomized data. This would need a more elaborate 
approach. 

6. Equal-weight edge randomization (EWER) 

Karsai et al. IPTD use a special null model that is de- 
signed for removing timing correlations between the contact 
sequences of adjacent edges, while retaining temporal char- 
acteristics associated with edges, including the distribution of 
inter-contact times on individual edges. Whole contact se- 
quences associated with edges, i.e. all contacts and their time 
stamps, are randomly exchanged between edges that have the 
same number of contacts. Thus single-edge patterns, such as 
burstiness, are retained, together with all properties retained 
by the RP model (number of contacts on each edge, network- 
level event frequency patterns, topological structure). This 
null model requires a large enough system so that there are 
enough edges with the same number of events. 

7. Edge randomization ( ER ) 

This null model is similar to the EWER model with the ex- 
ception that the sequences can be exchanged between edges 
that have any numbers of contacts. This corresponds to ran- 
domly exchanging the edge weights (measured as numbers of 



contacts) in the aggregated network spanned by the contact se- 
quence, and thus weight-topology correlations are destroyed 
in the null model. However, the inter-contact time distribu- 
tions of contact sequences of edges are not changed — the se- 
quences are just placed elsewhere in the network. 

8. Time reversal (TR) 

This null model is designed for assessing the frequency and 
importance of causal sequences J5) of contacts, where con- 
tacts trigger further contacts, and simply involves running the 
original event sequence backwards in time. If sequences of 
consecutive contacts would be caused by temporal correla- 
tions alone, similar numbers of such sequences should be ob- 
served when time runs forwards and backwards; the lack of 
such chains in the time-reversed null model compared to the 
original sequence can be attributed to "the arrow of time". 

9. Summary and guidelines 

The different randomized reference models discussed above 
retain and destroy specific kinds of topological and temporal 
correlations, and thus e.g. in studies of dynamical processes, 
they allow for pointing out the importance of various corre- 
lations: the most important correlations can be pinpointed by 
comparing the effects of different randomization models on 
the dynamics. The RE and RP models permute edges and 
contact times. Their simultaneous application (RE+RP) de- 
stroys all correlations except for patterns in the overall con- 
tact rate - this provides a good starting point for the limiting 
case of uncorrelated temporal networks, and by additionally 
randomizing contact times (RT), the overall patterns are also 
removed. When studying the roles of the exact contact tim- 
ings on edges and the correlations between adjacent edges, 
comparing the EWER and ER models to the RT model should 
do the trick, as the static network features are retained except 
for weight-topology correlations removed by the ER. 



VII. SPREADING DYNAMICS AND COMPARTMENTAL 
MODELS ON TEMPORAL GRAPHS 

One of the key insights of network theory is that the un- 
derlying network structure can strongly affect dynamic pro- 
cesses that are mediated by the edges. This is especially 
important for spreading processes, where biological or elec- 
tronic viruses, rumors, or pieces of information are transmit- 
ted through edges of physical contact, social ties, or electronic 
connections. For such processes, the network structure affects 
the speed of spreading as well as the extent of spread through 
the network through features such as short path lengths [ 161 1, 
the degree distribution IfTTl 11241 . degree correlations 1(171 . or 
community structure and correlations between tie strengths 
and network topology II 1 1 61 11221 . 

Many types of spreading processes can be modeled as com- 
partmental models, where each vertex is in one of the charac- 
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teristic states of the model II 1731 . i.e. belongs to one class or 
compartment ||3]|56). In the simplest SI model (Susceptible- 
Infective), individuals are initially susceptible, and become 
infected at some rate fj. when in contact with infective indi- 
viduals. In static networks, this process will eventually in- 
fect all vertices that can be reached from the source, and 
thus the network structure affects only the speed of spread- 
ing. The SIR model (Susceptible-Infective-Removed) shows 
richer dynamics. In this model, infected individuals recover 
and become immune to spreading, either with some probabil- 
ity /} per unit time or after a fixed period of time. The outcome 
of the SIR dynamics now depends on the interplay between 
the two rates (infection rate fi and recovery rate (3) as well as 
the network structure. The process may die out, only infect- 
ing local clusters of vertices, or percolate through the network 
in an epidemic fashion, such that a substantial fraction of in- 
dividuals is infected. Beyond these two simple models, more 
complicated versions have been formulated, such as the SIRS 
model, where immunity is not permanent but individuals may 
again become susceptible. 

Applying such models to static networks is equivalent to 
assuming that all interactions between vertices take place uni- 
formly in time. However, in reality this is usually not the case. 
The spreading of biological viruses depends on temporal pat- 
terns of contact, and the flow of information in social networks 
is influenced by the social activity patterns and rhythms of in- 
dividuals. Importantly, there is an increasing body of evidence 
of very large heterogeneity in the timings of such interactions, 
from burstiness in patterns of communication via electronic 
and physical mail 1171 l39l ETTl l69l 1 1 151 1 1571 . mobile telephone 
calls and text messages (21] [7T| 11071 11621 . instant messag- 
ing l95l . and patterns in proximity dynamics measured with 
RFID sensors [24, 641 11411 . Daily and circadian rhythms of 
human dynamics 1581 167ll94l 11021 [T03ll may also play a role. 
Furthermore, there are correlations where interaction events 
trigger further events, such as reception of an email triggering 
a forwarding event, or incoming calls causing outgoing calls 
in mobile call networks lETO 1711 11071 11201 . An overview of 
dynamic processes important for infectious disease spreading 
can be found in Bansal et al. J6). 

The above temporal inhomogeneities, together with the fact 
that spreading processes have to follow the time ordering of 
events, have important and at times drastic effects on the dy- 
namics of spreading on temporal graphs. In the recent years, 
such effects have been studied with the help of simulations 
building on empirical contact or interaction sequences, as well 
as with analytical tools. Although much remains unknown, 
some clear conclusions can already be drawn from the exist- 
ing literature, such as the importance of burstiness in slowing 
down epidemic-style spreading dynamics. 



A. Bursty event dynamics and slow spreading in 
communication networks 

Human communication dynamics is almost universally 
bursty, i.e. the timings between communication events devi- 
ate largely from uniform or Poissonian statistics and can in- 



stead be described with heavy-tailed or power-law distribu- 
tions. Consequently, the dynamics of spreading processes also 
differs from expectations. This failure of the Poissonian ap- 
proximation (that is typically assumed to hold for spreading 
models on static networks) was first addressed by Vazquez et 
al. II 1571 . who studied email activity patterns from university 
and service provider logs. 

Their study was motivated by the fact that the numbers of 
new infections by computer viruses and worms that spread 
through emails should decay exponentially, if the Poisson ap- 
proximation was accurate. However, reports of new infec- 
tions are typically still published years after the release of 
anti-viruses. In the Poisson approximation, the probability 
of one vertex to interact with another within a time interval 
dt is dtj{j), where (t) is the average interevent time. Thus, 
the time between consecutive contacts is exponentially dis- 
tributed, P(t) ~ exp(-f/(r)). However, for the email data 
studied by Vazquez et al., after correcting for finite observa- 
tion window, it was seen that P(f) ~ 1/t, followed by an 
exponential cutoff. 

The Poissonian and power-law inter-event time distribu- 
tions result in different generation times, that is, times be- 
tween one vertex becoming infected and transmitting the in- 
fection to its neighbor. On the basis of generation times, 
Vazquez et al. 015711 calculated an analytical approximation 
for the prevalence dynamics of the SI model assuming a tree- 
like structure, and showed together with simulations using 
empirical email sequences that the late-stage exponential de- 
cay of the number of new infections is significantly slower for 
the broad inter-event time distribution and matches well with 
data on real computer worms. 

Later, Min, Goh and Vazquez Ml 061 studied the SI model 
with power-law inter-event time distributions using theoret- 
ical arguments that assume a tree-like network structure, as 
well as with simulated spreading with uncorrelated power-law 
activity patterns and the priority-queue network model. They 
concluded that a power-law waiting time distribution leads to 
a power-law decay in the number of new infections in the long 
time limit, with an exponent determined through the genera- 
tion time distribution. 

More empirical evidence for slower-than-Poissonian 
spreading dynamics was provided by Iribarren and 
Moro ll6Tl l62l . who performed a viral marketing exper- 
iment with emails, where 31,183 individuals forwarded 
recommendations. They concluded that the large hetero- 
geneity found in the response times is responsible for the 
slow dynamics of information at the collective level. In 
their experiment, subscribers to an online newsletter were 
rewarded for recommending it via email to their friends. 
The viral spread of this recommendation email was tracked 
at every step, providing a detailed view on reception and 
forwarding events and their temporal correlations. As typical 
for such cascades, the average number of secondary cases per 
infected individual Rq was below the tipping point Rq = 1, 
in this case Rq « 0.26, and thus all cascades stopped in a 
finite number of steps. The cascades were observed to be 
treelike, with a very low clustering coefficient. It was seen 
that there is a large variability to the number of forwarded 
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recommendation emails and it was argued that this variability 
is not directly related to the degree of individuals in the 
email network. Furthermore, those individuals who became 
secondary spreaders, i.e. forwarded the recommendation, 
typically forwarded all their recommendation emails simul- 
taneously in a single spreading event, and did not remain 
spreaders for longer. The response times (times between 
reception and forwarding) were seen to follow a lognormal 
distribution. Traditional analytic epidemic models were seen 
to fail in predicting the speed and dynamics of the number 
of individuals who received and forwarded the message; 
according to the simple exponential growth equation, most 
new infections should happen during the early few days, 
whereas a significant fraction of new infections was observed 
even at the time scale of months. However, the observed 
slow dynamics was well captured by the non-Markovian 
Bellman-Harris branching model lf55l where the lognormal 
distribution of response times was used. 

Femandez-Gracia et al. [46 1 studied the effect of broad 
inter-event time distributions on the ordering dynamics of the 
Voter model; similarly to the slowing down of compartmen- 
tal spreading models, such distributions were observed to give 
rise to slow ordering dynamics of the model. 



B. Burstiness and other temporal and structural 
inhomogeneities 

Karsai et al. iPTD provided further insight into the effect 
of temporal heterogeneities on spreading dynamics by study- 
ing the behavior of the SI model. The model was simulated 
with real mobile telephone call and email contact event se- 
quences, together with reference models that destroy selected 



correlations (see Section VI Ci. The largest data source was 
the call database of a mobile telephone operator, containing 
about 325 million time-stamped call records over a period of 
120 days. For assessing the importance of different tempo- 
ral and structural inhomogeneities on the spreading dynam- 
ics, several reference models were used (with reference to the 
Section |VTC| >: RP, EWER, ER, a combination of RE and RP, 
and RT with independent Poisson processes on each edge. 
As in the above papers, inter-event times were seen to have 
a broad distribution; these were studied by calculating the dis- 
tributions for edges binned by event numbers and rescaling 
those by the average inter-event time in each bin, similarly to 
Candia et al. El and Miritello et al. II 1071 - It was also seen 
that the scaling breaks down for small times, around 20 sec- 
onds, indicating correlations and triggering of events. When 
compared to the original event sequence, the times to full 
prevalence (100% of vertices infected) were seen to be shorter 
for all null models. The RP model that destroys burstiness 
and correlations, except for daily patterns, gave faster spread- 
ing than the ER model that destroys weight-topology corre- 
lations but retains burstiness. Hence, the slowness of spread- 
ing can largely be attributed to the bursty event sequences on 
individual edges, in addition to the "weak-edge" bottleneck, 
i.e. edges between communities having a smaller contact fre- 
quency. However, the overall daily pattern where the system- 



level call frequency is low at night and peaks around noon and 
early evening was seen not to contribute significantly, when 
investigated with a Poissonian event-generating (RT) model. 
In addition, the EWER scheme that destroys triggering, i.e. 
temporal correlations between adjacent edges, was seen to re- 
sult in slightly slower-than-original spreading for short time 
scales, and slightly faster-than-original for long time scales. 

Interestingly, in contrast to these observations, Rocha et 
al. 013111 found that temporal order and correlations speed up 
epidemic spreading in their data set [ 1 30 1 of sexual contacts in 
Internet-mediated prostitution - the origins of this difference 
still remain unclear. Furthermore, in contrast to both findings, 
Stehle et al. found that when simulating the SEIR spreading 
dynamics on a temporal contact network of conference at- 
tendees recorded with proximity sensors of the SocioPatterns 
RFID platform [24|, the spreading dynamics is well described 
by a static aggregated network if the heterogeneity of the con- 
tact durations is taken into account as edge weights. 

Miritello et al. II 1071 also used mobile telephone call 
records (9 billion time-stamped calls of 20 million users over 
1 1 months) in their studies of simulated information spread- 
ing. They applied the SIR model with a homogeneous and de- 
terministic recovery time and variable transmission probabil- 
ity A with the empirical call sequence. The scaled distribution 
of relay times Tjj, that is, times between user i participating 
in a call with any other user and user i participating in a call 
with user j, was seen to be heavy-tailed while also display- 
ing a larger number of short relay times than expected from 
the Poissonian case. Similarly to Ref. [71], the abundance 
of short relay times was interpreted as a signature of group 
conversations, where calls trigger further calls and the activity 
patterns of adjacent edges are thus correlated. 

For SIR dynamics, such group conversations and corre- 
lated contact sequences were observed to give rise to larger 
spreading cascades for small values of A below the percola- 
tion point than for the time-shuffied reference model where 
the distribution of relay times approaches the Poisson distri- 
bution. However, for large values of A, the opposite behavior 
was observed. Thus, in the context of information propaga- 
tion, correlations from group conversations make information 
spreading more efficient at the local small scales when a more 
realistic low transmissibility A is used. Miritello et al. II 1071 
also proposed a way of mapping the dynamic SIR model to a 
static edge percolation model, similar to Newman 11 1 11 or 
Kenah and Robins l74l . They validated their approach by 
successfully predicting the percolation threshold for the SIR 
model on empirical data. This was done approximate the 
average number of secondary infections by an effective, dy- 
namic transmissibility — what they call the "dynamic strength 
of ties" — obtained from the mentioned mapping. 

It is worth noting that the speed of SI spreading — especially 
in the deterministic case where an infectious individual al- 
ways infects a susceptible individual upon contact — is related 
to temporal path lengths, latency and reachability ||59l 11201 . 
addressed elsewhere in this review. By placing constraints 
upon the timings between consecutive events defining a path, 
fastest temporal paths can also be representative of the path- 
ways taken by the SIR spreading process [ 120]. 
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Finally, in addition to temporal inhomogeneities slowing 
down spreading in human communication networks, a simi- 
lar effect has been observed for ants, however, compared to 
a different null model. In Ref. Ifl6l . temporal networks were 
constructed based on 30-minute video recordings of ants in 
6 colonies and tracking all contacts between individual ants. 
Compared to a kinetic null model, where ants were consid- 
ered as gas particles that randomly collide and change di- 
rections, Si-like information flow was observed to be signifi- 
cantly slower for the empirical contact sequences at long time 
scales. 

We also mention the modeling work by Kamp 11701 . where 
the author proposes a framework to study disease spreading 
in temporal networks, which simulates the contagion pro- 
cess without explicitly generating interval graphs. In Kamp's 
setup, vertices come in to the system with a degree sampled 
from a degree distribution (which then changes due to the evo- 
lution of the graph); they live for a limited time and their de- 
gree changes due to both the birth and death of neighbors, 
and rewiring of edges. This framework allows for some ap- 
proximate analytical treatments with generating functions and 
differential equation modeling but seems to require numerical 
simulations for a full characterization. 



C. Utilizing temporal structure for disease control 

The structure of contact patterns not only affects the spread- 
ing of disease, but this structure can also be exploited in con- 
trolling and preventing the spread. The most common pre- 
ventive intervention is vaccination that lowers the probability 
of people catching the disease and spreading it further. Typ- 
ically, one does not have to vaccinate the entire population 
to block the possibility of outbreaks. Already at some partial 
coverage / of vaccinees in the population, disease cannot any 
longer propagate. This effect is called herd immunity. Low- 
ering the threshold of herd immunity is an important goal in 
public health. In static network theory, there are methods that 
utilize the network topology to identify important targets for 
vaccination. Perhaps the most well-known is the neighbor- 
hood vaccination protocol of Cohen et al. [29 1 that works by 
repeating the following steps: 

1 . Take a random person in the population. 

2. Ask the person to name a friend (or rather someone that 
person meets regularly in such a way that the disease in 
question might spread). 

3. Vaccinate the friend. 

These steps are repeated until the desired fraction / of the 
population is vaccinated. The benefits of this scheme are 
twofold. First, it utilizes only local information — a person is 
expected to know his or her own contacts, not any third per- 
son's contacts — which really is a prerequisite rather than just 
an advantage. Second, it samples people in proportion to their 
degree k. The importance of a person with degree k with re- 
spect to disease spreading is, in classic network epidemiology, 
proportional to k 2 — roughly speaking, the probability that the 



person gets the disease is proportional to k, and the expected 
number of other people the person infects is also proportional 
to k (3). In a temporal network, it is important to remem- 
ber that future contacts cannot be used to determine the right 
targets for vaccination, and that past contacts cannot bene- 
fit from current vaccination, unless the contacts are repeated. 
The basic assumption of classic network epidemiology — that 
the network is static — means that the past contacts will always 
also persist into the future, and this does not hold in general 
in temporal networks. 

Lee et al. 1931 proposed an extension of the neighborhood 
vaccination scheme to temporal networks. They found that 
the strategy to ask about your most recent contact, or most 
frequent contact some time back in the past, improves the 
neighborhood vaccination in some real contact data sets (e.g. 
the prostitution sex network of Rocha et al. [ 1 30 1, the email 
data of Eckmann et al. 11391 . and the hospital proximity data 
of Liljeros et al. 11991 ). The response is different for differ- 
ent data sets, so for the hospital and prostitution data, the 
"most recent"-version is the most efficient, whereas for the 
email data the most frequent is more effective. This, Lee et 
al. argue, comes from the fact that the contacts along an edge 
in the hospital proximity and prostitution data are fairly lim- 
ited in time — two people who enter into a period of frequent 
contacts in either one of these data sets will rather likely be 
in contact with other persons a bit later. For the email data, 
the driving force behind the efficiency of the "most frequent"- 
protocol is that the contact frequency along an edge varies, 
but people in this dataset typically keep the relationship going 
throughout the data. From these simulations one can see that 
the temporal structure actually adds something that can be ex- 
ploited to local vaccination programs. We believe other types 
of population-sampling protocols that are affected by network 
structure in static simulations could need to be extended to 
temporal networks. 



VIII. FUTURE OUTLOOK 

In this review, we have illustrated how several systems can 
benefit from the temporal network approach, and discussed 
ways and methods of discovering and measuring network 
structure that lives in the time domain. Such structure be- 
comes important for network studies especially in the con- 
text of some dynamics taking place on the network: if there 
are inhomogeneities and correlations in the contact sequences 
between vertices, these inhomogeneities may have dramatic 
effects on dynamics mediated through the contacts. The con- 
ceptual differences between modeling dynamical processes on 
static and temporal networks deserve some attention. For ex- 
ample, for simple contact processes such as disease spread- 
ing models, in the static network approach one typically inte- 
grates two components in the model: when the contacts along 
an edge take place, and the probability for the disease be- 
ing transmitted during a contact between a susceptible and 
an infected vertex. The common assumption is that of con- 
tacts spread uniformly in time. In the temporal network ap- 
proach, the first component — the timings of the contacts — is 
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no longer a part of the spreading model, but rather an inte- 
gral part of the contact structure, i.e. the network itself. In 
general, for models of dynamical processes taking place on 
networks, moving from the static network framework to that 
of temporal networks is equivalent to removing the tempo- 
ral component related to contact timings from the model — 
such a component is part of almost every model, although it 
may not be explicitly visible — and considering temporal con- 
tact structure instead. But is the temporal network approach 
then nothing more than shifting a component of a dynami- 
cal model from the dynamical system to the underlying rep- 
resentation of interaction structure? What do we gain from 
such an approach, where the representation of the system in 
question necessarily becomes more complicated? First, as 
we have seen in this review, the temporal structure is crucial 
for studies of dynamical processes. Luckily, at least for pro- 
cesses involving humans, much of the data produced by to- 
day's technologies — be it mobile phone records 11711 11161 or 
wearable sensors 11241 l37l — comes in the form of contact se- 
quences and are thus directly suitable for temporal network 
studies. Thus one can go straight from the raw data to simu- 
lations of spreading processes (or some other dynamics), and 
analyze the role of the temporal and topological structure by 
comparing the results to reference models II 1 3 1 II - Second, 
one may learn a lot about the system and its driving forces 
by studying its temporal patterns and structure, from inter- 
contact statistics to patterns involving multiple vertices such 
as motifs. Additionally, temporal and topological structures 
can be correlated [114], and modeling this with an underlying 
static network structure is not straightforward. 

The study of temporal networks, their characteristic fea- 
tures, and their dynamics is still a rather young field, and there 
are many open questions and unexplored directions. Below, 
we list some of these issues: 

Generative models for temporal networks. There are only 
very few models for temporal networks and their contact se- 
quences, and one of the important open issues is clearly con- 
structing and studying parametrized, generative models of 
temporal networks, e.g. of human contact sequences with their 
characteristic features observed in real-world data, such as 
skewed inter-contact time distributions, bursty dynamics, and 
circadian and weekly rhythms. 

Measures for temporal network structure. Although a large 
number of measures and characteristics for temporal graphs 
have been discussed in this review, they have mostly been gen- 
eralizations of static network measures, and we feel that there 
is much room for improvements, e.g. in relation to simple 
measures of time-domain correlations of contact sequences. 
As an example, why does randomizing the order of contacts 
in a temporal network often, but not always, make spread- 
ing dynamics slower? It would be a little breakthrough if this 
question could be answered in terms of a simple, easily ob- 
servable measure (akin to the statement "high clustering co- 
efficient slows down disease spreading") II144I . In a sense, 
recent studies on temporal networks are going in the opposite 
direction from the early work on complex networks during the 
millennium's first years 111 121 — beginning with the effects of 
the temporal structure on dynamical systems, whatever this 



structure may be, rather than first quantifying and characteriz- 
ing the structural features of real systems, and only afterwards 
investigating the role of particular structural features on dy- 
namics. Today, we know that temporal network structure can 
make a difference, but not exactly how or why. 

Understanding the driving mechanisms. A third, largely 
unexplored theme is why contacts between two vertices 
in a temporal network happen when they happen — there 
are skewed inter-contact time distributions along individual 
edges, but in general, why? There are many papers about why 
two vertices become connected by an edge |112| and there 
are papers explaining time series of when a vertex does some- 
thing 0. However, presumably, what other vertices vertex i 
is connected to could affect the times when i does something, 
and vice versa. This question comes close to the goal of adap- 
tive network studies |48 | that model the feedback from net- 
work structure (and how it affects dynamics on the network) 
to the success of the agents forming the network (and how they 
seek to change their position in it). If one could include when 
contacts happen along an edge into adaptive network models 
and thereby explain some observed temporal-topological cor- 
relations, this would be a breakthrough (no matter what the 
objective system is). 

Inference problems. Another set of more statistics-related 
challenges concerns inference problems. How can one con- 
struct a temporal network from various amounts of informa- 
tion about the states of vertices or edges (Tj? How can one 
infer spreading chains, if one has an incomplete temporal net- 
work? Many inference problems on static networks should 
be rather different on temporal graphs as many fundamental 
properties — like the transitivity of edges or Menger's theo- 
rem 11511 — are only valid under rather strong assumptions. 

dynamical systems. If a system benefit from being modeled 
as a temporal network is not only a question about the struc- 
ture of the contacts, but also about the nature of the dynami- 
cal system acting over the contacts. Some dynamical systems 
might be more or less sensitive to temporal effects. In social 
information spreading, a person may, hypothetically, spread 
information only if he, or she, hears about it from two inde- 
pendent sources within a short period of time. Such a dynami- 
cal system should be sensitive to temporal effects like activity 
bursts and the order of events. 

Community, cluster or mesoscopic structure. The recent 
years have seen tremendous efforts for discovering meso- 
scopic structure in static networks in the form of communi- 
ties l42l . loosely defined as groups of vertices more densely 
connect within than between each other. Most of the litera- 
ture on community structure of static network focus on deriv- 
ing a method for decomposing the network from some kind 
of conceptually simple principle. Few, if any, studies seeks 
to identify structures known a priori to exist. The works in- 
corporating a time dimension into the community detection 
(like Refs. lfTT^[llll[T32l[T33ll ) operate on aggregated time- 
slices of the temporal network. One can imagine clustering 
algorithms based on more elaborate temporal structures, like 
time-respecting paths (a rare exception is Ref. [ 100|). 

Visualization. Visualization software is a great help for in- 
vestigating static networks. Even though small graphs are im- 
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possible to embed in Euclidean space such that the Euclidean 
distance between vertices is proportional to the graph distance 
(that would be ideal for a direct correspondence between the 
real graph — the terrain — and the visualization — the map), still 
one can make visualizations that capture much of the network 
structure. Typically they organize the vertices such that there 
is a positive correlation between the graph and Euclidean dis- 
tances |5Tj. This is enough to see differences between net- 
works of different degree-degree correlations 11 121 or to iden- 
tify dense clusters by the eye [42]. Stacking snapshots of a 
temporal network to a movie usually does a bad job to visu- 
alize temporal network structure, especially for sparse con- 
tact sequences. The time-line plots of e.g. Figs. [TJ [3] and 12 
resolve the vertices in one dimension, which makes it even 
harder to put vertices that are close in terms of latency (or 
some other distance-like metric for temporal networks) close 
to each other in Euclidean space. If one could find another 
clever way to layout a temporal network that captured reacha- 
bility and latency that would be a very valuable contribution. 

These open directions mentioned above mainly concern 
theoretical and methodological developments. However, the 
real acid test of temporal networks as a fruitful paradigm is its 
application to concrete, specific problems in population biol- 
ogy, cell biology, ecology, neuroscience, social and political 
sciences, economics, chemistry and so forth. So far, the tem- 



poral network framework has, with not so many exceptions, 
been investigated theoretically rather than used to explain the 
world around us. Yet most complex systems in the world are 
time-dependent, dynamical and in motion. As we believe that 
the temporal networks framework is really a tool for advanc- 
ing science, we hope to see theoretically-minded researchers 
bringing it into collaborations with their applied colleagues. 
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